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ABSTRACT 

We compute and investigate four types of imprint of a stochastic background of primordial magnetic fields (PMFs) on the cosmic microwave 
background (CMB) anisotropies: the impact of PMFs on the CMB temperature and polarization spectra, related to their contribution to cosmo¬ 
logical perturbations; the effect on CMB polarization induced by Faraday rotation; the impact of PMFs on the ionization history; magnetically- 
induced non-Gaussianities and related non-zero bispectra; and the magnetically-induced breaking of statistical isotropy. We present constraints 
on the amplitude of PMFs derived from different Planck data products, depending on the specific effect that is analysed. Overall, Planck data 
constrain the amplitude of PMFs to less than a few nanogauss, with different bounds depending on the considered model. In particular, indi¬ 
vidual limits coming from the analysis of the CMB angular power spectra, using the Planck likelihood, are Mpc < 4.4 nG (where 6, Mpc is 
the comoving field amplitude at a scale of 1 Mpc) at 95 % confidence level, assuming zero helicity. By considering the Planck likelihood, based 
only on parity-even angular power spectra, we obtain Sim P c < 5.6 nG for a maximally helical field. For nearly scale-invariant PMFs we obtain 
Bi Mpc <2.0 nG and B\ m pc < 0.9 nG if the impact of PMFs on the ionization history of the Universe is included in the analysis. From the analysis 
of magnetically-induced non-Gaussianity we obtain three different values, corresponding to three applied methods, all below 5 nG. The constraint 
from the magnetically-induced passive-tensor bispectrum is BiMpc < 2.8 nG. A search for preferred directions in the magnetically-induced passive 
bispectrum yields B t Mpc < 4.5 nG, whereas the the compensated-scalar bispectrum gives B t Mpc < 3 nG. The analysis of the Faraday rotation of 
CMB polarization by PMFs uses the Planck power spectra in EE and BB at 70 GHz and gives B i Mpc < 1380 nG. In our final analysis, we consider 
the harmonic-space correlations produced by Alfven waves, finding no significant evidence for the presence of these waves. Together, these results 
comprise a comprehensive set of constraints on possible PMFs with Planck data. 
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1. Introduction 

1.1. Cosmic magnetism 

Magnetic fields are one of the fundamental and ubiquitous com¬ 
ponents of our Universe. They are a common feature of many 
astrophysical objects, starting from the smallest up to the largest 
observed scales (for reviews see Ryu et al. 2012 and Widrow 
et al. 2012). In particular, large-scale magnetic fields are ob¬ 
served in almost every galaxy, starting from the Milky Way, with 
possible hints of their presence also in high-redshift galaxies 
(Beck, 2000; Bernet et al., 2008; Wolfe et al., 2008), suggest¬ 
ing an early origin for the galactic fields. 

Large-scale magnetic fields are also probed in galaxy clus¬ 
ters, both through the measurement of the Faraday rotation effect 
on the light of background galaxies and through radio emission 
from the halos and relics of the clusters (Govoni & Feretti, 2004; 
Ferrari et al., 2008; Feretti et al., 2012). These large-scale mag¬ 
netic fields have measured amplitudes that range from a few to 
several microgauss. Recent Faraday rotation measurements from 
low-density intercluster regions suggest the presence of large- 
scale magnetic fields also in cosmic structure filaments (Neronov 
et al., 2013a). A recent addition to the large-scale magnetic field 
constraints comes from the interpretation of Fermi-LAT data 
(Neronov & Vovk, 2010). High-energy y-rays (in the TeV band) 
emitted by blazars generate electron-positron pairs when inter¬ 
acting with the optical and IR background light. These pairs re¬ 
emit at lower energies (GeV), but Fermi -LAT data do not show 
this re-emitted flux in the GeV band. One possible explanation 
for this observation is a deflection of the electron-positron pairs 
due to the presence of a diffuse magnetic field (see Neronov 
& Vovk 2010, Taylor et al. 2011, Tavecchio et al. 2011, Vovk 
et al. 2012, and Neronov et al. 2013b for details on this expla¬ 
nation and see Broderick et al. 2012 for alternative scenarios). 
Constraints on the GeV emission provide lower limits on the 
amplitude of intergalactic fields of the order of 10 1 s -10 15 G 
(Tavecchio et al., 2010; Taylor et al., 2011; Dermer et al., 2011; 
Vovk et al., 2012), if this scenario is correct. 

The origin of large-scale magnetic fields is strongly debated. 
Several mechanisms have been proposed and one popular hy¬ 
pothesis is that the observed large-scale fields are remnants of 
fields that existed from the earliest times, i.e., primordial fields. 
During structure formation the adiabatic compression and tur¬ 
bulent shock flows would naturally lead to an amplification of 
initial seeds (which may act in addition to astrophysical mech¬ 
anisms of large-scale magnetic field generation, like AGN ejec¬ 
tion and galactic dynamos; for reviews see Widrow 2002 and 
Giovannini 2004b). A Kolmogorov-like magnetic power spec¬ 
trum has been observed in the central region of the Hydra cluster 
(Kuchar & EnBlin, 2011), supporting the idea that the observed 
extragalactic magnetic fields are largely shaped and amplified by 
hydrodynamic al processes. Primordial magnetic fields (PMFs) 
can naturally provide the initial seeds to be amplified into the ob¬ 
served large-scale fields. Several early-Universe scenarios pre¬ 
dict the generation of cosmological magnetic fields, either dur¬ 
ing inflation (Ratra, 1992), with a suitable breaking of conformal 
invariance of electromagnetism (Turner & Widrow, 1988), dur¬ 
ing phase transitions motivated by particle physics (Vachaspati, 
1991; Grasso & Riotto, 1998), or via other physical processes 
(Durrer & Caprini, 2003; Ichiki et al., 2006). The importance 
of PMF studies lies not only in the possibility of PMFs being 
the progenitors of the observed cosmic magnetic fields, but also 
in them providing a new potential observational window to the 
early Universe (for reviews see Kahniashvili 2005, Giovannini 
2008, Kunze 2013, and Durrer & Neronov 2013). 


1.2. Imprints of primordial magnetism 

PMFs leave imprints on several cosmological observables and 
can be constrained with different cosmological data sets. In par¬ 
ticular, interesting constraints come from their influence on Big 
Bang nucleosynthesis, which provides upper limits of the or¬ 
der of 0.1 fiG (Grasso & Rubinstein, 1995; Kahniashvili et al., 
2010), and from their impact on large-scale structure formation 
(see, e.g., Shaw & Lewis 2012 and Fedeli & Moscardini 2012). 

Another limit on PMFs based on Big Bang Nucleosynthesis 
(BBN) is derived by Caprini & Durrer (2002). Gravitational 
waves can be produced by PMFs, in particular before neutrino 
free streaming. The upper limit on the amount of gravitational 
waves allowed at nucleosynthesis to not spoil the BBN predic¬ 
tions poses a constraint on the amplitude of PMFs, which is es¬ 
pecially strong for causal magnetogenesis mechanisms. 

It has been suggested that PMFs could have an influence on 
the formation of the filamentary large-scale structure, in the pi¬ 
oneering work of Wasserman (1978) and by Kim et al. (1996). 
In particular, Battaner et al. (1997) conclude that magnetic fields 
with comoving strengths lower than 1 nG have negligible effect 
and that for strengths higher than lOnG the influence is too high 
to be compatible with observations. 

Other constraints on PMFs come from their impact on the 
thermal spectrum of the CMB radiation. PMFs can induce spec¬ 
tral distortions of both early and late type via injection of dissi¬ 
pated magnetic energy into the plasma through damping pro¬ 
cesses (Jedamzik et al., 2000; Sethi & Subramanian, 2005; 
Kunze & Komatsu, 2014). PMFs imply also photon emission 
and absorption through the cyclotron process, which, depend¬ 
ing on the amplitude of PMFs, could in principle play a role in 
the generation and evolution of spectral distortions and in par¬ 
ticular in the thermalization process (for a general overview see 
Burigana et al. 1991, Hu & Silk 1993, and Chluba & Sunyaev 
2012). On the other hand, the cyclotron process is relevant only 
at very long wavelengths. Thus, for realistic shapes of distorted 
spectra and PMFs with amplitudes compatible with current con¬ 
straints, the cyclotron contribution is found to be much less 
important than radiative Compton and bremsstrahlung contri¬ 
butions in re-establishing a blackbody spectrum (Burigana & 
Zizzo, 2006) and also extremely small in generating a polarized 
signal (Zizzo & Burigana, 2005). Thus, the current constraints 
derived from COBE-FIRAS present good prospects for possible 
future observations (ledamzik et al., 2000; Kunze & Komatsu, 
2014; Chluba et al., 2015) while polarization anisotropies di¬ 
rectly induced by PMFs are not significantly affected by the cy¬ 
clotron process associated with PMFs. 

Many of the stronger and more robust constraints come from 
the impact of PMFs on the CMB anisotropies. CMB data are a 
crucial source of information for investigating and constraining 
PMF characteristics, understanding their origin, and exploring 
the possibility of them being the seeds that generated the ob¬ 
served large-scale magnetic fields. 

The first release of Planck 1 data in 2013 has led to 
some of most stringent constraints on PMFs (see Planck 
Collaboration XVI, 2014). The scope of this paper is to pro¬ 
vide the “Planck constraints on PMFs” through combined anal- 


1 Planck (http://www.esa.int/Planck) is a project of the 
European Space Agency (ESA) with instruments provided by two sci¬ 
entific consortia funded by ESA member states and led by Principal 
Investigators from France and Italy, telescope reflectors provided 
through a collaboration between ESA and a scientific consortium led 
and funded by Denmark, and additional contributions from NASA 
(USA). 
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yses of the temperature and polarization data. The results of 
this paper are derived from Planck products, which are based 
on the work done in Planck Collaboration I (2015), Planck 
Collaboration II (2015), Planck Collaboration III (2015), Planck 
Collaboration IV (2015), Planck Collaboration V (2015), Planck 
Collaboration VI (2015), Planck Collaboration VII (2015), and 
Planck Collaboration VIII (2015). Before going into the descrip¬ 
tion of our analysis, we briefly discuss the most important PMF 
models. 

The simplest PMF model is that of a homogeneous field. 
This model cannot be included in a homogeneously and isotrop¬ 
ically expanding cosmological model. It needs to be anal¬ 
ysed in the context of an anisotropic cosmological model with 
associated isotropy-breaking predictions (Kahniashvili et ah, 
2008), and has already been strongly constrained by COBE data 
(Barrow et ah, 1997). More recently, Adamek et ah (2011) have 
reconsidered the impact of a homogeneous large-scale magnetic 
field on the CMB anisotropies, with the addition of the contribu¬ 
tion from free-streaming particles like neutrinos. The presence 
of the anisotropic neutrino stress induces a compensation of the 
magnetic field effect, allowing for a larger homogeneous-field 
amplitude than concluded from previous analyses. 

The most widely used model of PMFs is a stochastic 
background modelled as a fully inhomogeneous component 
of the cosmological plasma (we neglect PMF energy density 
and anisotropic stress contributions at the homogeneous level), 
with the energy momentum tensor components (quadratic in 
the fields) on the same footing as cosmological perturbations. 
Within this model, PMFs leave several signatures in the CMB 
temperature and polarization anisotropy patterns, inducing also 
non-Gaussianities. 

In this paper, we predict and analyse four different types 
of PMF signatures: the impact on the CMB power spectra in 
temperature and polarization; the impact on polarization power 
spectra induced by Faraday rotation; the impact on CMB non- 
Gaussianities and the related non-zero magnetically-induced 
bispectra; and finally the impact on the statistics of CMB 
anisotropies and in particular the breaking of statistical isotropy, 
given by induced correlations in harmonic space. 

The energy momentum tensor of PMFs sources all types of 
cosmological perturbations, i.e., scalar, vector, and tensor per¬ 
turbations. Magnetically-induced pertubations have some cru¬ 
cial differences with respect to the primary perturbations. Firstly, 
PMFs generate vector perturbations that are sourced by the 
Lorentz force and, unlike the primary ones, are not decay¬ 
ing. Secondly, magnetically-induced perturbations are not sup¬ 
pressed by Silk damping (Hu & White, 1997; Subramanian & 
Barrow, 1998a). Their impact on the CMB temperature power 
spectrum is dominant on small angular scales, where the primary 
CMB is suppressed, and therefore they can be strongly con¬ 
strained by high-resolution CMB data. Of additional interest are 
helical PMFs, which may be generated during inflation through 
mechanisms like pseudoscalar coupling and may be crucial for 
connecting these PMFs to the magnetic fields observed on large 
scales. A helical PMF generates parity-violating correlations, 
such as TB and EB (Caprini et ah, 2004; Kahniashvili & Ratra, 
2005; Kahniashvili et ah, 2014), which can be used to constrain 
PMFs with CMB polarization data. After recombination, PMFs 
are dissipated via two additional effects that take place in the 
magnetized, not fully ionized, plasma: ambipolar diffusion and 
magnetohydrodynamical (MHD) turbulence. The dissipation in¬ 
jects magnetic energy into the plasma, heating it and thereby 
modifying the optical depth of recombination, with an impact on 


the primary CMB power spectra (Sethi & Subramanian, 2005; 
Kunze & Komatsu, 2014, 2015; Chluba et ah, 2015). 

PMFs have another effect on the primary polarization 
anisotropies. They induce Faraday rotation, which rotates E- 
modes into /1-modes, thus generating a new /l-mode signal, and 
vice versa. This signal grows with decreasing observational fre¬ 
quency and therefore is a good target for Planck's low-frequency 
channels. 

PMFs modelled as a stochastic background have a non- 
Gaussian contribution to CMB anisotropies, even if the mag¬ 
netic fields themselves are Gaussian distributed, since the com¬ 
ponents of the energy momentum tensor are quadratic in the 
fields and therefore approximately follow^ 2 statistics (Brown & 
Crittenden, 2005). In particular, PMFs generate non-zero higher- 
order statistical moments. The third-order moment, the CMB 
bispectrum, can be used as a probe to derive constraints on 
PMFs that are complementary to the previously mentioned ones. 
Planck polarization data thus provide a new way of probing 
PMFs, namely through the magnetically-induced polarization 
bispectrum. 

The presence of PMFs induces and sustains the propaga¬ 
tion of Alfven waves. These waves have an impact on the 
statistics of the CMB anisotropies and in particular induce spe¬ 
cific correlations between harmonic modes (Kahniashvili et ah, 
2008). It is possible to use this effect to constrain the ampli¬ 
tude of Alfven waves and thus indirectly constrain the PMF 
amplitude (Durrer et al., 1998; Kim & Naselsky, 2009; Planck 
Collaboration XXIV, 2014). 

In the MHD limit, assuming that the fields are only modified 
by cosmic expansion, the magnetic field strength decreases as 
a~ 2 , where a is the cosmological scale factor. Throughout, we 
will use a “comoving” magnetic field, defined as B = a 1 B (phys) , 
where B <phys > is the physical strength of the magnetic field. 

1.3. Structure of the paper 

The paper is structured as follows. In Sect. 2 we describe the 
analysis of the impact of helical and non-helical PMFs on CMB 
power spectra in temperature and polarization and we derive the 
constraints on the PMF amplitude and spectral index coming 
from Planck data. In Sect. 3 we present three different analyses 
of the magnetically-induced bispectrum, specifically two analy¬ 
ses of the magnetically-induced passive bispectrum and an ana¬ 
lytical treatment of the magnetically-induced scalar bispectrum. 
In all cases we derive the constraints on the amplitude of PMFs 
with a scale-invariant spectrum using Planck non-Gaussianity 
measurements. In Sect. 4 we present our analysis of the Faraday 
rotation signal induced by PMFs and we derive constraints from 
Planck low frequency polarization data. In Sect. 5 we present 
the analysis of the impact of Alven waves on statistical correla¬ 
tions in harmonic space and the associated constraints on Alfven 
waves derived from Planck data. We summarize our conclusions 
in Sect. 6 


2. Impact of primordial magnetic fields on the 
CMB power spectra 

PMFs affect cosmological perturbations and may leave signifi¬ 
cant imprints on the CMB power spectra in temperature and po¬ 
larization. Accurate prediction of these signatures allows us to 
derive constraints on PMF characteristics from CMB anisotropy 
data from Planck using the Planck likelihood. In this section we 
derive the predictions for the magnetically-induced power spec- 
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tra in temperature and polarization, considering helical and non¬ 
helical PMFs, and present the resulting constraints on PMFs. 

2.1. Magnetic modes 

When considering a stochastic background of PMFs, we can ne¬ 
glect the contribution of energy density and anisotropic stress at 
the homogeneous level. The magnetic energy momentum ten¬ 
sor can be seen as describing perturbations carrying energy den¬ 
sity and anisotropic stress and inducing a Lorentz force on the 
charged particles of the plasma. PMFs source all types of pertur¬ 
bations; scalar, vector, and tensor. In the past years, several dif¬ 
ferent analyses of magnetically-induced perturbations have been 
performed. Some examples from the wide literature of the field 
concern magnetically-induced scalar perturbations (Giovannini, 
2004a; Kahniashvili & Ratra, 2007; Yamazaki et al., 2007, 2008; 
Finelli et ah, 2008; Giovannini & Kunze, 2008c, a,b; Bonvin 
& Caprini, 2010; Bonvin, 2010; Kunze, 2011), while other 
treatments also include magnetically-induced vector and ten¬ 
sor perturbations (Subramanian & Barrow, 1998b; Durrer et ah, 
2000; Kahniashvili et ah, 2001; Mack et ah, 2002; Caprini & 
Durrer, 2002; Subramanian & Barrow, 2002; Subramanian et ah, 
2003; Lewis, 2004; Caprini, 2006; Paoletti et ah, 2009; Shaw 
& Lewis, 2010). We can identify three different classes of intial 
conditions for magnetically-induced perturbations; compensated 
(Giovannini, 2004a; Finelli et ah, 2008), passive (Lewis, 2004; 
Shaw & Lewis, 2010), and inflationary (Bonvin et ah, 2013). 

In this paper we focus on the two magnetically-induced 
modes that are present for all types of PMFs produced prior 
to decoupling, independent of their generation mechanism, i.e., 
the compensated and passive modes. We do not consider spe¬ 
cific inflationary initial conditions (Bonvin et ah, 2013) to main¬ 
tain the generality of the PMFs we constrain. For the same rea¬ 
son we neither consider a possible cross-correlation between the 
magnetically-induced and the adiabatic mode motivated by in¬ 
flation (Jain & Sloth, 2012). 

2.1.1. Compensated modes 

The compensated modes are the regular magnetically-induced 
modes. These are the regular (finite at r —> 0) solutions of the 
perturbed Einstein-Boltzmann equations, including the magnetic 
contributions after neutrino decoupling. These modes are called 
“compensated” because the magnetic contributions to the metric 
perturbations in the initial conditions are compensated by fluid 
modes to leading order. The initial conditions are the solutions of 
the Einstein-Boltzmann equation system for large wavelengths 
at early times, with the perturbed quantities expanded in power 
series of kr (where k is the perturbation wavenumber and r is 
the conformal time). When performing this calculation for the 
magnetically-induced modes, the growing regular mode requires 
the source terms in the equations for the metric perturbations to 
vanish at the lowest order. This can only be realized by a com¬ 
pensation between the magnetic terms and the perturbed quanti¬ 
ties of the fluid. 

2.1.2. Passive modes 

The second class, the passive modes, is generated by the pres¬ 
ence of a PMF before neutrino decoupling. Without neutrinos 
free-streaming, there is no counterpart in the fluid to balance 
the anisotropic stress of the PMF. This generates a logarithmi¬ 
cally growing mode (in conformal time, which diverges for early 


times). After neutrino decoupling, the anisotropic neutrino stress 
compensates the anisotropic stress due to the PMF, leading back 
to the compensated case described before. But an imprint of this 
logarithmically growing mode survives neutrino decoupling in 
form of a constant offset on the amplitude of the inflationary 
non-magnetic mode (the primary cosmological pertubations of 
the standard model without PMFs). This amplitude offset is due 
to the continuity condition for the matching of the initial condi¬ 
tions before and after neutrino decoupling. Passive modes have 
a logarithmic dependence on the ratio between the neutrino de¬ 
coupling time and the generation time of the PMF, i.e., their 
amplitudes grow as h(k) oc In (t v /t b ) (where r v is the neutrino 
decoupling time and t b is the PMF generation time). The pas¬ 
sive modes, unlike the compensated ones, evolve following the 
standard non-magnetic equations and influence only scalar and 
tensor perturbations. 

2.2. Impact of non-helical PMFs on the CMB angular 
power spectra 

Our analysis is based on previous treatments of magnetically- 
induced compensated and passive scalar, vector, and tensor 
modes presented by Lewis (2004), Finelli et al. (2008), Paoletti 
et al. (2009), and Shaw & Lewis (2010). 

At linear order, PMFs evolve like a stiff source and we can 
therefore discard the back-reaction of gravity onto the stochas¬ 
tic background of PMFs. Prior to the recombination epoch the 
electric conductivity of the primordial plasma is very large. We 
therefore consider the limit of infinite conductivity, in which the 
induced electric field is zero. In this limit, the temporal evolution 
of the PMF reduces to f? (phys) (x, r) = B(x)la{r) 2 , where B(x) is 
the comoving field. 2 

We model a non-helical stochastic PMF with a power-law 
power spectrum, where the two-point correlation function is de¬ 
scribed by 3 

(Bfk) B)(k')) = ^ 5 & (k - k') (< 5ij - hkj) P B (k ), (1) 

where P B (k) = A B k" B and k, denotes a cartesian component of 
a normalized wave vector. In this model, the PMF is character¬ 
ized by two quantities, the amplitude of the power spectrum, A B , 
and the spectral index n B . The latter is one of the major discrim¬ 
inating factors between generation mechanisms, since different 
mechanisms generate fields with different spectral indices (for 
example, causal mechanisms generate fields with n B > 2; Durrer 
& Caprini 2003). While magnetically-induced compensated per¬ 
turbations do not suffer from Silk damping, PMFs are neverthe¬ 
less suppressed on small scales by radiation viscosity (Jedamzik 
et al., 1998; Kahniashvili et al., 2001). To account for this damp¬ 
ing we introduce a sharp cut-off in the PMF power spectrum at 
the damping scale k d. 


2 We choose the standard convention in which the scale factor is 
a(r 0 ) = 1 at the present time r 0 . 

3 For the Fourier transform and its inverse, we use 

Y(k,r ) = J d 3 xe' k x Y(x,r ), 

m ’ T) = f H e ^m,r), 

where Y is a generic function. 
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For the amplitude we use the convention to smooth over a 
comoving scale of A — 1 Mpc, 





dkk 2 

2lT 2 


-k 2 ,Y 


P B (k ) = 


Ab 

4 7 r 2 / (n s + 3 


r 


n B + 3 \ 


( 2 ) 


For the damping scale we use (Subramanian & Barrow, 1998a; 
Mack et al„ 2002) 4 


k D = (5.5 x lO 4 )^ 


j(n. h h 2 \" 


<sr( 


2 n 

U/Mpcj 


n B + 3 




h "b+ 5 


0.022 j 


1 , 1=1 Mpc 


Mpc~', 


(3) 


where the I1 (X) are the scalar, vector, and tensor components of 
the energy momentum tensor, P tJ = d, ; - k,kj. and we sum over 
repeated indices. Such convolutions can be written in terms of 
spectra as 

(n* iS \k)n (S \k')) = |n (S, (k)| 2 s(k-kf), 

= l -\u (W \k)\ 2 Pij(k)6(k-k'), 

(nf\k)n™(k'j) = l -\TF\k)\ 2 M m (k)6(k-k r ), 

where Ad, ;; ; = P„P t \ + PaPjt - PijPti- With this convention, the 
relevant components of the energy momentum tensor become 


where h is the reduced Hubble constant, Hq = 
100/t kins 1 Mpc -1 , and Ob is the baryon density parameter. 5 
Magnetically-induced scalar, vector, and tensor perturbations 
are sourced by the energy momentum tensor components due to 
PMFs, together with the Lorentz force contribution. The energy 
momentum tensor of the PMFs is 


B 2 (x) 

87ra 4 (r) ’ 

(4) 


(5) 

y - b,(x) /rVxj), 

(6) 


where the components are all quadratic in the magnetic field. 
The power spectra of the perturbations are therefore fourth-order 
in the magnetic field and given by convolutions of the magnetic 
power spectrum. The two-point correlation function of the spa¬ 
tial part of the energy momentum tensor is 6 

( K a b (k) K cd (k')) = J' d< ^J S ab Scd 


I p B (k)\ 2 = — f d 3 pP B (p)P B (\k-p\)(l+p 2 ), (8) 

10247T 3 J n 

= T2sk?£ d>pP ‘ (p) 

xP B (\k - p\) [l + p 2 + 4 jP(jP - p)\ , (9) 

\n <V) (k)\ 2 = -±— s ^d 3 pP B (p)P B (\k-p\) 

xf(l+/? 2 )(l - r 2 ) + y/^- r /?)] , (10) 

| n<T) W| 2 = ^ \/ P P B ip)P B (\ k -p\) 

x(l+2 y 2 +rp 2 ), (11) 

where p - p ■ (k - p)/\k - p\, y - k ■ p. [5 = k ■ (k - p)/\k - p\, 

and f l denotes the volume with p < k n . 

The conservation equations for the fields give a relation be¬ 
tween the scalar projection of the anisotropic stress, the energy 
density, and the Lorentz force, cr B = y + L B 1 (which reduces 
to a simple relation between Lorentz force and anisotropic stress 


(B,(q) B,(k - q) B m (p) B m (k' - p)) 


f 


d 3 q d 3 p 
327T 5 


(■ B a (q ) B h (k - q) B c (p) B d (k’ - pfj . 

We can then obtain scalar, vector, and tensor correlation func¬ 
tions. 


(n * (S) (k) n (S1 (k')) = 6 ab 6 cd (,d ab (k ) Kcd (k ')), 

(n ; (Vi (k) n ( p(k')) = k a p ih(k) k' c p Jd (k') (,«* ab (k) K cd (k ')), 


(n;j T> (k) n' ; T) (A:')) = 


p ia (k) p jb (k) - -p u (k)p ab (k) 


p u: (k')P ld (k')- -p tl (k')P cd (k') 


{< lh (k)K cd (k')), (7) 


4 Note that, unlike in previous analyses (Finelli et al., 2008; Paoletti 
et al., 2009; Paoletti & Finelli, 2011), we explicitly include the depen¬ 
dence of the damping scale on the baryon density to account for the 
updated cosmology with respect to the previous treatments. 

5 For a quasi-scale-invariant PMF power spectrum ( n B ~ -3), the es¬ 
timated value for k D varies slightly, by about 30-40 %, for different ap¬ 
proaches (Subramanian & Barrow, 1998a; Jedamzik et al., 2000; Kunze 

& Komatsu, 2015). However, this does not affect the main results of this 
paper significantly. 

6 We use the convention that Latin indices run from 1 to 3, while 
Greek indices run from 0 to 3. 


for vector modes, If^ 1 = k'L* V) ). This relation simplifies the 
treatment, reducing the number of correlators to be computed by 
a factor of 2. We use the analytic solutions to the convolutions 
derived by Finelli et al. (2008) and Paoletti et al. (2009) for fixed 
spectral indices. For general n B , we use the fits to the generic an¬ 
alytic solutions provided by Paoletti & Finelli (2011). These fits 
simplify the computation due to the presence of hypergeomet¬ 
ric functions in the analytic solutions. The infrared behaviour of 
the spectra, which is relevant for CMB anisotropies, depends on 
the spectral index. In particular, the spectra describe white noise 
for indices greater than n B = — whereas they are infrared- 
dominated, as k 2 " s+3 , for smaller indices. We use the initial con¬ 
ditions derived by Lewis (2004), Paoletti et al. (2009), Paoletti & 
Finelli (2011), and Finelli & Paoletti (2015) for scalar and com¬ 
pensated tensor modes, vector modes, and passive modes, re¬ 
spectively. We use an extended version of the CAMB code (Lewis 
& Challinor, 2011) that includes all magnetic contributions to 
calculate predictions for the CMB power spectra in temperature 
and polarization. 

2.2.1. Compensated modes 

In Fig. 1, we show the predictions for magnetically-induced 
compensated modes. This shows that the dominant compen¬ 
sated contributions to the angular power spectra are given by 


7 Note that we use the notation of Ma & Bertschinger (1995) for the 
scalar anisotropic stress a B . 
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Fig. 1. Magnetically-induced CMB TT (top left), TE (top right), EE (bottom left), and BB (bottom right) power spectra. The solid 
lines represent primary CMB anisotropies, the dotted lines represent magnetically-induced compensated scalar modes (except for 
the BB panel, where it represents the lensing contributions and the solid line represents primary tensor modes with a tensor-to-scalar 
ratio of r = 0.1), the dashed lines represent vector modes, whereas the dot-dashed lines represent magnetically-induced compensated 
tensor modes. We consider PMFs with fiiMpc = 4.5 nG and = -1. 


, 3 , 

CM 



Fig. 2. Dependence of the magnetically-induced CMB power spectrum on the spectral index. For all plotted cases, the amplitude 
is B i Mpc = 4.5 nG. The black lines show primary CMB anisotropies; for the other colours we refer to the legend. Left: scalar 
contributions, right : vector contributions. 


the scalar and vector modes. In particular, because magnetically- 
induced perturbations are not suppressed by Silk damping, a sig¬ 


nificant contribution of magnetically-induced modes arises on 
small angular scales, where the primary CMB fluctuations are 
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suppressed. As we will show, the impact of PMFs on the CMB 
power spectrum at high multipoles is particularly relevant for 
the high-precision Planck data, allowing us to derive strong con¬ 
straints on the PMF amplitude. Since magnetically-induced per¬ 
turbations are solely sourced by energy momentum tensor com¬ 
ponents due to PMFs, the shape of the magnetically-induced 
spectra strongly depends on the PMF spectral index. In Fig. 2, 
we show this dependence for the temperature power spectrum 
for scalar and vector perturbations. We note the qualitatively dif¬ 
ferent dependence for two regimes. For n B > -3/2 the angular 
power spectrum remains flat (i.e., 1(1 + 1 )Cf oc f 2 ) with a rescal¬ 
ing of the amplitude due to the amplitude of the Fourier spectra, 
whereas for ng < -3/2 the shape varies according to the infrared 
domination of the energy momentum tensor of the PMFs. 


2.2.2. Passive modes 

In addition to compensated initial conditions, we also consider 
passive tensor modes. The magnetically-induced passive modes 
are not completely determined by the amplitude and spectral in¬ 
dex of the PMF, but depend also on the ratio t v /t b . For tensors 
we specifically have li(k) oc n (T) (k) ln(r y /r B ), where h{k) is the 
tensor metric perturbation. This ratio may vary between 10 17 and 
10 6 for fields originating at the grand unification energy scale 
(GUT) and at later phase transitions, (Shaw & Lewis, 2010). We 
consider the passive tensor modes, since, for red spectra, they 
give the dominant contribution on large angular scales, where 
the compensated modes are subdominant. In Fig. 3, we com¬ 
pare the magnetically-induced passive tensor modes for a nearly 
scale-invariant spectrum, n B = -2.9, for PMFs generated at the 
GUT scale (for which t v /t b = 10 17 ) with the corresponding 
dominant compensated modes. In Fig. 4, we show the depen¬ 
dence of the CMB power spectrum due to passive tensor modes 
on the spectral index, for GUT scale PMFs, and its dependence 
on the time ratio, showing the two extreme values of the possible 
range. We note that the compensated vector modes dominate at 
small angular scales, but the passive tensor modes give a con¬ 
tribution at low and intermediate multipoles for red spectra. For 
bluer spectra, the passive spectrum becomes steeper and there¬ 
fore subdominant with respect to primary CMB fluctuations on 
large angular scales and with respect to vector modes on small 
angular scales. A similar behaviour can be observed in the scalar 
passive mode, which has the same origin as the tensor one but 
in the scalar sector. Shaw & Lewis (2010) have shown that, just 
like the tensor mode, the scalar pasive mode becomes relevant 
on large angular scales for nearly scale-invariant power spectra. 
We will show below that the dominant passive contribution to 
the constraints on the PMF amplitude is given by tensor modes. 


2.3. Impact of helical PMFs on CMB anisotropies 

In addition to the ubiquitous presence of magnetic fields in the 
Universe, astrophysical observations show that some galaxies 
might have a helical magnetic field structure (Widrow, 2002; 
Vallee, 2004). Following the hypothesis that the amplification 
of PMFs may have played a role in the generation of large-scale 
magnetic fields, the observed magnetic helicity may be related 
to helicity of the PMFs. A helical intergalactic magnetic field 
could also be related to a possible CP violation recently hypoth¬ 
esized by Tashiro et al. (2013) in an indirect study of cosmo¬ 
logical large-scale magnetic fields in voids using secondary y- 
ray data. Helicity of the PMFs influences MHD processes in 
the early plasma, as well as cosmological perturbation dynam¬ 


ics, allowing different processes of energy transport, for exam¬ 
ple the inverse cascade mechanism (Biskamp, 2003). These pro¬ 
cesses of energy transport play a role in the early time evolu¬ 
tion of the PMFs and may have an impact on our understanding 
of their generation mechanisms, especially if PMFs with a non¬ 
zero helicity are generated. Moreover, the presence of helicity 
would test possible modifications of Maxwell’s theory by con¬ 
straining parameters describing the gauge invariance (i.e., mass 
of the photon) and Lorentz invariance (i.e., existence of a pre¬ 
ferred frame of reference) as discussed by Carroll et al. (1990). 
Thus it would carry information about particle physics at very 
high temperatures (above 1 TeV). 

A possible way to detect magnetic helicity directly from 
CMB data is to study the polarized CMB (cross-) power 
spectra. A non-zero helicity in the PMFs changes the ampli¬ 
tudes of the parity-even power spectra and induces parity-odd 
cross-correlations between the temperature and B-polarization 
anisotropies and E- and B-polarization anisotropies (Pogosian 
et al., 2002; Caprini et al., 2004; Kahniashvili & Ratra, 2005; 
Ballardini et al., 2015). Such parity-odd cross-correlators are 
also generated by Faraday rotation, but only due to homoge¬ 
neous PMFs, not due to the stochastic model for PMFs consid¬ 
ered in this paper (Kosowsky et al., 2005). Parity-odd signals 
may therefore give a more direct possibility for studying helical 
PMFs. 

Our present study for a stochastic background of helical 
PMFs is an extension of the model described in the previous sec¬ 
tion. We consider the impact of such PMFs on CMB anisotropies 
in temperature and polarization. Using the exact expressions for 
the energy-momentum tensor components including the helical 
contribution discussed by Ballardini et al. (2015), we derive the 
predictions for the impact of a helical PMF on the CMB power 
spectra in temperature and polarization. 

The most general ansatz for the two-point correlation func¬ 
tion, built on Eq. (1), but taking into account an antisymmetric 
part (Pogosian et al., 2002), is 

(B,(k) B)(k')) = ( ^6 0) (k-k') 

\(S,j ~ kjkj) P B (k) + i etji h P H (k)\ , (12) 

with Pfi(k) = A H k"" . From a geometrical point of view P B (k) 
denotes the symmetric part and Pn(k ) the antisymmetric part of 
the correlator. The totally antisymmetric tensor is related to 
the parity violation under a transformation k —» -k. The realiz¬ 
ability condition gives \Pn{k)\ < P B (k). 

In this case the model is described by four parameters, two 
amplitudes and two spectral indices, where A B is the ampli¬ 
tude of the power spectrum of the fields, the same as defined 
in Eq. (2), and Ah is the amplitude of the power spectrum of the 
helical part of the PMFs. These amplitudes can be expressed in 
terms of mean-square values of the magnetic field and of the he¬ 
lical component, respectively. As done for the amplitude of the 
magnetic field in Eq. (2), we can express the amplitude of the 
helical component on a comoving scale of /l as (Ballardini et al., 
2015) 


si 



dkk 3 
In 1 


e 


-k 2 P 


\P H (k)\ = 


\Ah\ p 
4n 2 A' lH+i 


n H + 4 1 


■ (13) 


The helical spectral index needs to satisfy tin > -4 for conver¬ 
gence. 

The helical term in Eq. (12) generates new sources that 
contribute to the energy momentum tensor. All components of 
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Fig. 3. Magnetically-induced CMB TT (top left), TE (top right), EE (bottom left), and BB (bottom right) power spectra due to 
passive tensor modes, compared with the ones due to compensated modes. The solid lines represent primary CMB anisotropies, 
the dotted lines represent magnetically-induced compensated scalar modes (except for the BB panel, where it represents the lensing 
contribution), the dashed lines represent vector modes, whereas dot-dashed lines represent magnetically-induced passive tensor 
modes. We consider PMFs with B\ \ 1pc = 4.5 nG and n B = -2.9. 
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Fig. 4. Dependence of the magnetically-induced CMB power spectrum due to passive tensor modes on the spectral index for a GUT- 
scale PMF (left) and comparison between the two extremes for the time ratio t v /tb (right). The black lines show the primary CMB 
anisotropies; for the other colours we refer to the legend. Solid lines represent PMFs generated at the GUT scale, t v /tb = 10 17 , 
whereas dashed lines represent PMFs generated at late times, t v /t b = 10 6 . 
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Fig. 5. CMB TT (top left ), TE (top right), EE (bottom left), and BB (bottom right) power spectra due to helical PMFs compared 
to the ones due to non-helical PMFs. Solid lines are non-helical predictions, while dashed lines are helical predictions. Blue are the 
scalar modes, green the vector, and red the compensated tensor modes. We consider PMFs with B\ Mpc = 4.5 nG and n K = -1. 


the energy momentum tensor are quadratic in the fields, so the 
Fourier-space two-point correlation function generates symmet¬ 
ric and antisymmetric sources: symmetric ones due to products 
of Pb with Pb (the components found in the non-helical case) 
and products of P H with P H ; antisymmetric ones, which gener¬ 
ate odd-parity angular power spectra, from products of Pg with 
P H . The Fourier components of the energy momentum tensor 
are 


I Ps(k)\ z 


l/ J /4^')| n0 n helical 

1 


512 


A 


dpP H (p) P H (\k-p\)p, 


(14) 


\Lf(k)\ 2 = |Z®(*)| 


1 


non-helical 
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64 7t 2 a 8 


x d 3 p p P H (p) P H (\k -p\)(p- 2yjS) , (15) 


K’wf = |n‘ v> (*)L_ h „ 

1 


512 


tX‘ 


cP p P„(p) P H (\k-p\) (p-yP) , (16) 


\n m (k)\ 2 = |tf T, (*)| 


non-helical 

■3. 


128 7T 5 I‘ 


d pP H (p) Prt(\k - p\)y(J, 


(17) 


where the non-helical parts are given by Eqs. (8)— (11). The sys¬ 
tem of equations for the background and the perturbations, as 
well as the initial conditions, are unmodified. 

In Fig. 5 we show the predictions for the magnetically- 
induced compensated modes, considering the additional contri¬ 
butions to the sources due to helicity. We consider the maximally 
helical case. Ah = Ag, with equal spectral indices, n H = n H , us¬ 
ing the solutions to the helical energy momentum tensor compo¬ 
nents derived by Ballardini et al. (2015). The predictions show 
no difference in the shape and in the slope of the angular power 


spectra, with a small shift in the amplitude, which is always 
smaller for the helical case, at least for scales relevant for the 
CMB. Magnetic helicity induces parity-odd cross-correlations 
between the E- and B-polarization anisotropies, as well as be¬ 
tween temperature and B-polarization anisotropies. The parity- 
odd cross-correlations are sourced by the mixed terms in the 
correlation function of the energy momentum tensor, propor¬ 
tional to f d 3 p Pg(p) Ph(\Ic - p\). These terms, after decompo¬ 
sition, contribute to the vector and tensor sources as 

k ( v ) ® i 2 = \/p{ p ^p)PH(\k- P \) 

[/?(i -Y L )-(yP-p)y\ 

+P'h(p) Ps(\k - p\) [y(l -p 2 ) - (yp - p)p \J ,(18) 
\A m (k)\ 2 = -±- 5 ^d 2 p\p B (p)P H {\k-p\)[li{\+y 1 )\ 

+P H (p)PB(\k-p\)[y(l+f3 2 )}y (19) 

In the limit of small momenta, i.e., k <sc k n , the spectra never 
show a white-noise behaviour, |/F X) (k)| = const, contrary to what 
happens for the symmetric and non-helical parts. For ng + n H > 
-2 the behaviour is proportional to kk'^ +r,H+2 and for ns + < 

-2 it does not depend on the damping scale and is proportional 
to k" B+nH+2 . 


2.4. Constraints from the CMB temperature and 
polarization power spectra 

Here we present the constraints from Planck on helical and non¬ 
helical PMFs. In the literature there are several previous stud¬ 
ies that already derived constraints on PMFs using different 


9 













































Planck Collaboration: Constraints on primordial magnetic fields 


combinations of observed CMB power spectra (Caprini, 2010; 
Yamazaki et al., 2010; Paoletti & Finelli, 2011; Shaw & Lewis, 
2012; Paoletti & Finelli, 2013; POLARBEAR Collaboration 
et al., 2015). We use an extended version of the CosmoMC code 
(Lewis & Bridle, 2011), modified to include the magnetic con¬ 
tributions to the CMB power spectra, as described in the previ¬ 
ous subsections, and to include the parameters characterizing the 
PMFs in the Markov chain Monte Carlo analysis. 

We assume a flat Universe and a CMB temperature Tq = 
2.7255 K, and we use the BBN consistency condition (Ichikawa 
& Takahashi, 2006; Hamann et al., 2008). We restrict our anal¬ 
ysis to three massless neutrinos. A non-vanishing neutrino mass 
would not modify the results since it would only enhance the 
power on large scales in the presence of PMFs for the compen¬ 
sated modes, where the PMF contribution is less relevant (Shaw 
& Lewis, 2010). The pivot scale of the primordial scalar is set to 
k,, = 0.05 Mpc 1 . We consider the lensing effect for the primary 
CMB power spectrum and follow the method implemented in 
the Planck likelihood to marginalize over astrophysical residuals 
and secondary anisotropy contamination of the small-angular- 
scale data (Planck Collaboration XI, 2015). This contamination 
is particularly relevant for the PMF scenario, since PMFs im¬ 
pact mainly small angular scales. If this contamination is not 
properly considered it may lead to biased constraints on PMFs 
(Paoletti & Finelli, 2013). We sample the posterior using the 
Metropolis-Hastings algorithm (Hastings, 1970), generating be¬ 
tween four and sixteen parallel chains and imposing a conser¬ 
vative Gelman-Rubin convergence criterion (Gelman & Rubin, 
1992) of/? - 1 < 0.01. 8 We vary the baryon density cj b = Qb h 2 , 
the cold dark matter density a> c = £l c h 2 , the reionisation optical 
depth T re i on , the ratio of the sound horizon to the angular diame¬ 
ter distance at decoupling 6, the scalar amplitude ln(A s 10 10 ), and 
the scalar slope « s . In the MCMC analysis, we include the mag¬ 
netic parameters B\ Mpc and n B for the compensated modes, and 
add the parameter T rat = t v /t b whenever we also consider the 
passive tensor mode. We use flat priors for the magnetic param¬ 
eters in the ranges [0,10] for/?! Mpc /nG, [-2.9,3] for n B (n B > -3 
to avoid infrared divergence in the PMF energy momentum ten¬ 
sor correlations). We sample r rat logarithmically, with a flat prior 
on log 10 r rat in the range [4,17], 

2.4.1. Likelihood 

We derive the constraints on PMFs using the Planck likelihood, 
which is described in detail in Planck Collaboration XI (2015). 
Here we give a brief summary of the main points. The Planck 
likelihood is based on the Planck 2015 data and considers both 
temperature and polarization. As in 2013, we use a hybrid ap¬ 
proach with the combination of two likelihoods, one dedicated 
to low / and the other to high /. 

The Planck low-/ likelihood is a fully pixel-based likelihood 
with temperature and polarization treated jointly and at the same 
resolution, A' s j ( | c = 16. The Grange is 2 < / < 29 in TT, TE, 
EE, and /?/?. The likelihood is based on the foreground-cleaned 
LFI maps at 70 GHz and the temperature map derived by the 
component separation method Commander using 94% of the 

x This convergence criterion is based on the analysis of the variance 
within a chain and between chains. The posterior marginal variance is 
a weighted average of the different variances. If all the chains have 
reached convergence, this value will be very close to the variance within 
each single chain. The parameter R is defined as the square root of the 
ratio of the posterior marginal variance and the variance within a chain. 
With this definition, the closer R is to unity, the closer to convergence 
the chains are. 


sky at frequencies from 30 to 353 GHz (Planck Collaboration 
IX, 2015). The polarization map covers 54% of the sky and 
is derived from the 70 GHz Q and U maps cleaned with the 
30 GHz map as a synchrotron template and the 353 GHz map as 
a dust template (see Planck Collaboration XI, 2015). This like¬ 
lihood is denoted as “lowP” throughout the paper. Contrary to 
the 2013 analysis, where a combination of Planck temperature 
and WMAP9 polarization data was used, the 2015 low-/ likeli¬ 
hood is based entirely on Planck data for both temperature and 
polarization. 

The Planck high-/ likelihood is based on a Gaussian 
approximation (Planck Collaboration XV (2014) and 
Planck Collaboration XI (2015) for polarization) and cov¬ 
ers the /-range 30 < / < 2500. It uses the half-mission 
cross-power spectra of the 100 GHz, 143 GHz, and 217 GHz 
channels, measured in the cleanest region of the sky far from 
the Galactic plane and bright point sources. The sky fractions 
considered are 66 % of the sky for 100 GHz, 57 % for 143 GHz, 
and 47% for 217 GHz in temperature, whereas in polarization 
they are 70%, 50%, and 41%, respectively. The likelihood 
takes foregrounds and secondary anisotropies into account. In 
particular, for the temperature spectra it considers the contri¬ 
butions of dust, clustered Cosmic Infrared Background (CIB), 
thermal and kinetic Sunyaev Zeldovich effect (tSZ and kSZ), 
the cross-correlation between tSZ and CIB, and a Poissonian 
term for unresolved point sources for the temperature spectra. 
In polarization, only the dust contribution is considered. Each 
model is parameterized as a template contribution to the Q with 
a free amplitude. The dominant contribution for / < 500 is dust, 
whereas high /-modes are dominated by point sources and in 
particular the CIB for the 217 GHz auto-correlation. For the 
details of the foreground modelling see Planck Collaboration 
XI (2015). This high-/ likelihood is denoted as “Planck TT", 
for temperature only, or “Planck TT, TE, EE”, for temperature 
plus polarization, throughout the paper. 

2.4.2. Constraints with compensated scalar and vector 
contributions 

We perform an analysis with the Planck 2015 baseline likeli¬ 
hood. In Table 1 we report the derived constraints. The constraint 
on the PMF amplitude is Z?i m P c < 4.4 nG at the 95 % confidence 
limit (CL) for the case that includes temperature and polarization 
data both at low and high multipoles. The same constraint results 
when including the polarization only at low /. As in previous 
analyses, PMFs with positive spectral indices are constrained to 
lower amplitudes than PMFs with negative spectral indices. In 
Fig. 6 we present the results of this analysis compared with the 
Planck 2013 constraints (Planck Collaboration XVI, 2014). The 
upper limits from current Planck data are slightly higher than 
those obtained in 2013 (Planck Collaboration XVI, 2014), where 
the constraint from Planck data alone was /?im P c <4.1 nG. The 
weaker constraint can be explained by several changes of the 
2015 data and likelihood with respect to 2013. The change in 
the calibration (see Planck Collaboration I 2015), the different 
likelihood implementation, and the different models for the fore¬ 
ground residuals are all factors that contribute to the changed 
upper limit. In fact, all these factors, including also the slightly 
higher spectral index n s with respect to 2013, are mimicking a 
slightly larger signal in the temperature anisotropies, which is 
compatible with larger values of the PMF amplitude. 

We now include the polarization data in the analysis. 
Although the impact of PMFs on TE and EE polarization is less 
important than on temperature anisotropies, we show the results 
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Table 1. Mean parameter values and bounds of the central 68 % 
CL from Planck TT,TE,EE (left column) and Planck TT (right 
column). When consistent with zero, the upper bound of the 
95 % CL is reported. Note that // (l is a derived parameter. The 
posterior of the spectral index ub is strongly prior-dependent 
since B\ vi pc is consistent with zero. 


Parameter Planck TT,TE,EE + lowP Planck TT + lowP 


cj b . 0.0222 ± 0.0002 0.0222 ± 0.0002 

cj c . 0.1198 ±0.0015 0.1197 ±0.0022 

6 . 1.0408 ± 0.0003 1.0408 ± 0.0005 

r reion . 0.078 ±0.017 0.075 ±0.019 

log[A s 10- 9 ]_ 3.09 ±0.03 3.08 ±0.04 

n s . 0.963 ±0.005 0.964 ± 0.007 

Ho . 67.77 + _ 0 0 ™ 67.82 + _™* 

^iMpc/ n G. < 4.4 < 4.4 

n B . <-0.008 <-0.31 


Table 2. Upper bounds of the central 95 % CL for the PMF am¬ 
plitude. C stands for compensated mode, C+P for compensated 
plus passive modes, T re ion prior indicates the case where instead 
of the low-f polarization likelihood, as a cross-check, we used a 
Gaussian prior on the optical depth, r re i 0n = 0.07 + 0.02. 


6iMpc/nG 


TT,TE,EE+\owP\ C . <4.4 

7T+lowP: C. <4.4 

TT,TE,EE+lowP: C+P ... <4.5 

7T+lowP:C+P. <4.5 

TT + r reion prior: C+P. < 4.4 


for the case of Planck data, which include also high-C TE and 
EE polarization, in Fig. 6. Although the shape of the posterior 
changes slightly, there is no net improvement on the 95 % CL 
upper bound on B\ vipc with the addition of the high-7 TE and 
EE polarization. 

2.4.3. Constraints with passive tensor contributions 

As described in the previous subsection, in addition to regu¬ 
lar magnetically-induced compensated modes, the presence of 
PMFs prior to neutrino decoupling generates passive modes. 
In Fig. 3 we show that the passive tensor modes may give 
the dominant magnetic contribution to the CMB power spec¬ 
tra for a nearly scale-invariant PMF power spectrum. Their in¬ 
clusion in the analysis may therefore be relevant for the con¬ 
straints on PMFs. We include the passive tensor contribution in 
the MCMC code with the addition of the parameter r mt , with 
the settings described above. We perform MCMC analyses with 
the Planck 2015 likelihood, combining the low-7 temperature 
and polarization data either with high-7 temperature data or with 
high-7 temperature and polarization data, i.e., 7’7 +lowP and 
7T,7’£',£ , £ , +lowP. 

Figure 7 and Table 2 present the results, compared to the re¬ 
sults of the case that includes only compensated contributions. 
The 95% CL constraint on the PMF amplitude is B\ Mpc < 
4.5 nG, which implies that the addition of the passive tensor con¬ 
tribution does not improve the constraint on the amplitude of the 
PMFs. This result is expected on the basis of the shape of the an¬ 
gular power spectra due to passive tensor modes and their strong 
dependence on the PMF spectral index. This mode is basically 




n B 


Fig. 6. Comparison of the constraints on the smoothed PMF am¬ 
plitude (top) and the spectral index ( bottom ) from the 2015 tem¬ 
perature and temperature plus polarization data with the 2013 
results for magnetically-induced compensated initial conditions 
only. 


a primary tensor mode with and amplitude that depends on the 
PMFs. Its spectrum flattens on large angular scales and then de¬ 
cays on intermediate ones for red PMF spectra, whereas it ac¬ 
quires a steeper shape for blue PMF spectra, but with a much 
lower amplitude than for compensated vector modes. Therefore, 
passive tensor modes only contribute significantly for nearly 
scale invariant indices. In Fig. 8 we present the two-dimensional 
plot for the PMF amplitude and the spectral index. It shows the 
strong degeneracy between the two parameters, meaning that the 
same magnetically-induced power spectrum can be realized with 
different pairs of amplitude and spectral index. The effect of pas¬ 
sive tensor modes on the CMB temperature angular power spec- 
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tram is dominant over the primary CMB anisotropies only for 
a very limited range of spectral indices, near the scale-invariant 
case. Therefore, the degeneracy between amplitude and spectral 
index reduces the influence of the passive tensor mode contribu¬ 
tion on the constraints on the amplitude and we do not see any 
improvement in adding the passive tensor mode. While the con¬ 
straint on the amplitude is almost unchanged, the spectral index 
in the MCMC analysis is sensitive to the contribution of the pas¬ 
sive mode. In the lower panel of Fig. 7 we show the different 
shapes of the posterior distributions for the PMF spectral index 
for different data combinations. The inclusion of the passive ten¬ 
sor mode influences the low spectral index part of the posterior, 
while the compensated modes influence the high spectral index 
part. In addition to the analyses that consider the combination of 
passive and compensated modes, we perform an analysis includ¬ 
ing only the contribution of the passive tensor mode. We obtain 
B i Mpc < 6.5 nG at 95% CL using Planck TT+ lowP. This re¬ 
sult shows that the contribution of the passive term alone, when 
considering the spectral index and the generation epoch as free 
parameters, does not have the constraining power of the combi¬ 
nation of passive and compensated modes. 



2.4.4. Impact of astrophysical residuals 

We adopt the Planck likelihood treatment of astrophysical 
contaminants as described in the Planck likelihood paper. 
Considering the complexity of this model and the number 
of nuisance parameters involved, we investigate whether fore¬ 
ground residuals have an impact on the constraints on PMFs. 
Specifically, we investigate possible degeneracies with the fore¬ 
ground parameters by considering the two-dimensional distri¬ 
butions of the magnetic parameters and foreground parameters. 
The relevant cases are shown in Fig. 9. There, we present the 
two-dimensional distributions of the PMF amplitude and the 
Poissonian amplitudes for the three frequencies considered in 
the Planck high-f likelihood: 100 GHz, 143 GHz, 217 GHz, and 
the 143 x217 GHz cross-spectrum. We note that especially for 
the 143 GHz and the 143x217 GHz analyses, a weak degener¬ 
acy between the two parameters is seen. This result may indi¬ 
cate an impact of the astrophysical residual modelling on the 
PMF constraints. To investigate this issue, we perform an analy¬ 
sis by fixing the four parameters for the Poissonian amplitudes to 
their best-fit values from the standard Lambda Cold Dark Matter 
model. This analysis yields a limit of B\ yi pc < 3.0 nG, which 
is smaller than the constraint obtained in the case where the pa¬ 
rameters associated to astrophysical residuals are free to vary. 
This result has no statistical significance, but demonstrates that 
there is an impact of the astrophysical residuals on the PMF con¬ 
straints when the data considered require a complex model for 
the residuals. The shape of the dominant PMF contributions to 
the angular power spectrum on small angular scales is responsi¬ 
ble for this degeneracy. In fact, the steep slope of the vector mode 
may be degenerate with astrophysical residual contributions, as 
shown by Paoletti & Finelli (2013). For comparison we show an 
analogous plot for the 2013 analysis in Fig. A.l, which consid¬ 
ers the degeneracy for the foreground parameters of the Planck 
2013 likelihood. We note how, in contrast to the 2015 analysis, 
there is only a small degeneracy with the Poissonian amplitude 
at 143 GHz. This result shows the importance of the foreground 
residual modelling for the PMF constraints. 

In contrast to the Poissonian terms, we do not observe any 
degeneracy with the other foreground components, as shown in 
Fig. A. 2, including the clustering component of the foreground 
residuals, which is the other dominant contribution on small an- 



n B 

Fig. 7. Constraints on the smoothed PMF amplitude (top) and 
spectral index (bottom) from Planck temperature data with and 
without the passive tensor contribution. Constraints including 
both compensated and passive modes are indicated with C+P 
in the legend, constraints using only compensated modes are 
marked with C. 

gular scales at the frequencies considered in this analysis. The 
fact that we do not observe a degeneracy in this case is due to 
the difference in the spectral shape of the PMF contribution and 
the clustering term. Although both are relevant on small angular 
scales the slightly different shapes break the degeneracy. 

2.4.5. Constraints for specific PMF models 

Planck 2015 results confirm what has been observed in previous 
analyses, namely that CMB data allow negative PMF spectral 
indices with larger field amplitudes than positive indices. The 
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0 2 4 6 


^1 Mpc [llC] 

Fig. 8. PMF amplitude versus the spectral index for the baseline 
Planck 2015 case. C+P denotes the case where both compen¬ 
sated and passive modes are considered, whereas C indicates the 
case with only compensated modes.The two contours represent 
the 68 % and 95 % confidence levels. 
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Fig. 9. Two-dimensional posterior distributions of the PMF am¬ 
plitude versus the parameter describing the Poissonian term of 
unresolved point sources for the three frequencies considered in 
the likelihood. The two contours represent the 68 % and 95 % 
confidence levels. 


spectral index of the PMFs is the main discriminating factor 
among possible generation mechanisms. 

Some cases are of particular interest due to their connec¬ 
tion with specific classes of generation mechanisms. In particu¬ 
lar, PMFs generated during phase transitions or via second order 
perturbative effects, vector perturbations, etc., are characterized 
by positive spectral indices, equal or greater than 2. To inves¬ 
tigate the maximal amplitude allowed by Planck data for fields 
of this type, we perform two dedicated analyses, the first with 
fixed index n B = 2, and the second only restricted to positive 
spectral indices for PMFs. We include both compensated and 
passive modes in the analysis, giving B\ m pc < 0.011 nG at 95 % 
CL (B \ Mpc < 0.012 nG at 95% CL, when considering Planck 


Table 3. 95 % CL upper bounds of the PMF amplitude for fixed 
spectral index with compensated plus passive tensor modes. 


n B . 

2 10-1 

-1.5 

-2 

-2.5 

-2.9 

^lMpc/flG • • 

. 0.011 0.1 0.5 3.2 

4.8 

4.5 

2.4 

2.0 


TT ,T E ,EE+\oviP) for the n B = 2 case and Bim P c < 0.55 nG at 
95 % CL for n B > 0. 

We consider a third case of interest, the almost scale- 
invariant fields with n B = -2.9. This specific case is connected to 
PMF generation from inflation and is studied to test the strength 
of the passive tensor modes in constraining the amplitude of 
the PMFs. Moreover, we want to compare the results obtained 
from the Planck power spectra with those coming from the non- 
Gaussianity analysis, presented in the next section, which is per¬ 
formed for this spectral index as well. We obtain B\ m pc <2.0 nG 
at the 95 % CL. Note that this nearly scale-invariant case is dom¬ 
inated by the tensor passive mode. In fact, when we consider 
only the tensor passive contributions, excluding the compensated 
ones, we obtain the same result as in the passive-compensated 
combined case, B\ m pc < 2.0 nG at the 95 % CL. 

Together with the passive tensor mode there is also a scalar 
passive mode, as shown by Shaw & Lewis (2010). When we 
include this scalar passive contribution in our analysis we obtain 
again B\ m P c < 2.0 nG at the 95 % CL. We can therefore conclude 
that the passive scalar contribution is subdominant with respect 
to the tensor one. As will become clear in the next section, this 
result shows that the constraining power of the angular power 
spectrum is comparable to the one of the non-Gaussianity. 

In addition to these specific types of PMFs we have per¬ 
formed some analyses with fixed spectral index, choosing a grid 
of values covering the full range we sample. With respect to the 
case where the spectral index is a free variable, the cases with 
fixed spectral index are expected to give stronger constraints on 
the PMF amplitude, with a trend in agreement with the general 
case, because one of the two parameters describing the PMF is 
fixed. In Table 3 we present the results of these analyses. We note 
how, as expected, the trend of the results with fixed spectral in¬ 
dex is in agreement with the one of the two-dimensional plot of 
Fig. 8 obtained with the generic sampling of the index. We note 
also how the constraint is weakest for n B = -1.5 as expected 
from the impact on the angular power spectrum. 

2.4.6. Constraints from the BICEP2/ Keck-P/anck joint 
analysis 

We perform an analysis using the recent BICEP2/ Keck-Planck 
cross-correlation (indicated as BKP; BICEP2/Keck Array and 
Planck Collaborations, 2015) in addition to the Planck 2015 
data. The BKP likelihood is obtained from the BB and EE band- 
powers for all cross-spectra between the BICEP2 /Keck maps and 
the Planck maps at all frequencies (BICEP2/Keck Array and 
Planck Collaborations, 2015). We consider both compensated 
and passive contributions and study two cases, one in which we 
leave the spectral index free to vary and another in which we 
fix it to n B = -2.9. In the latter case, there is a contribution to 
the fi-mode polarization on large angular scales from the passive 
tensor mode. Figure 10 shows the comparison of the results of 
these two analyses with the results obtained from Planck data 
alone. The constraints are B ] \ 1pc < 4.7 nG for the case with free 
spectral index and Bim P c < 2.2 nG at the 95 % CL for the case 
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-61 Mpc [^G] 


Fig. 10. Probability distributions for the PMF amplitude includ¬ 
ing the BICEP2I Keck-Planck cross-correlation, compared with 
the one based only on Planck data. Top : the case in which the 
spectral index is free to vary, bottom : the case with n B — -2.9. 


with n b = -2.9. These are slightly higher upper bounds for the 
amplitude of PMFs, but they are fully compatible with the re¬ 
sults derived from Planck data alone. We note that the posterior 
distribution for the nearly scale-invariant case changes with the 
addition of the BKP data but it does not show any significant 
deviation from the posterior based only on Planck data. 

2.4.7. Constraints with maximally helical contributions 

We perform an MCMC analysis including the maximally heli¬ 
cal contribution. We restrict our analysis to the case of tempera¬ 
ture and polarization with only even cross-correlations. The odd 
cross-correlations TB and EB are present only in the lowP like- 



Fig.ll. PMF amplitude constraint for the helical case (solid 
black) compared with the non-helical case (dashed red). The dot¬ 
ted blue line shows the constraint on the amplitude of the helical 
component as an alternative interpretation of the constraints on 
the amplitude of PMFs with a helical component. 


lihood, therefore only for very low multipoles where the signal 
from helical PMFs is negligible. Thus we do not include odd 
cross-correlators in our analysis. 

We perform an analysis using the PlanckTT-v lowP likeli¬ 
hood. The constraint on the PMF amplitude in the maximally 
helical case is B ] \ lpc < 5.6 nG at the 95% CL. A comparison 
with the corresponding results for the non-helical case is shown 
in Fig. 11. The analysis with the Planck TT ,TE,EE+\owP like¬ 
lihood gives B i Mpc < 5.8 nG at the 95% CL. As in the non¬ 
helical case, the inclusion of high-f polarization does not im¬ 
prove the constraints. Figure 5 shows that magnetic fields with 
a maximally helical component produce smaller CMB fluctua¬ 
tions in temperature and polarization than non-helical fields of 
the same strength. As a result of this, the amplitude of max¬ 
imally helical magnetic fields are less constrained than non¬ 
helical fields for this Planck 2015 data release. When consid¬ 
ering helical PMFs, we have two components that contribute to 
the magnetically-induced perturbations, as shown in Eq. (12), 
a symmetric and an antisymmetric part, represented by P B and 
Ph , respectively. These power spectra can be associated with 
two amplitudes of the field, B\ vi pc associated with the symmet¬ 
ric part and Si Mpc associated with the antisymmetric part (see 
Eq. 13). In the maximally helical case the two amplitudes are 
not independent from each other, they are related through the 
conditions A H = A B and n B = n H . Therefore we constrain a sin¬ 
gle amplitude, which can be expressed either through fiiMpc or 
Si Mpc - The constraint B\ M pc <5.6 nG can thus be converted into 
the constraint Si m P c < 4.6 nG at the 95 % CL. Figure 11 shows 
the posterior distribution for the amplitude expressed as Si Mpc 
in blue. 
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2.4.8. Constraints from the impact of PMFs on the CMB 
anisotropies via their impact on the thermal history 
of the Universe 

Primordial magnetic fields are damped on scales smaller than the 
photon diffusion and free-streaming scale. This leads to heating 
of ordinary matter (electrons and baryons), which affects both 
the thermal and ionization history of the Universe (Subramanian 
& Barrow, 1998a; Jedamzik et al., 2000; Sethi & Subramanian, 
2005; Schleicher et ah, 2008; Kunze & Komatsu, 2014; Chluba 
et ah, 2015), leading to a Compton-y distortion of the CMB and 
changes in the CMB power spectra through modifications of the 
Thomson visibility function around decoupling. 

Two heating mechanisms have been discussed in the liter¬ 
ature, one due to decaying magnetic turbulence at very small 
scales and the other due to ambipolar diffusion (e.g., Sethi & 
Subramanian, 2005). In this paper, we follow the approach de¬ 
scribed by Chluba et ah (2015) to incorporate these heating 
mechanisms. 9 

We perform an analysis considering the combination of 
the heating terms with the gravitational contribution of PMFs. 
Considering ambipolar diffusion, decaying magnetic turbulence, 
and gravitational effects we obtain an upper limit of B\ M pc < 
0.90 nG at 95 % CL for nearly scale-invariant PMFs with n B = 
-2.9. We obtain the same result, namely B\ yi pc < 0.90 nG at 
95 % CL, when dropping the gravitational effect and consider¬ 
ing only the impact of PMFs on the primary CMB anisotropies 
through their heating effect. These results show that the dom¬ 
inant contribution is given by the heating terms. We have also 
performed analyses with the two terms of ambipolar diffusion 
and decaying magnetic turbulence considered separately. The 
results show that the two terms are roughly at the same level 
in constraining PMFs, with a slightly stronger contribution from 
the decaying magnetic turbulence term (see also Chluba et al., 
2015). 

Together with the Planck TT + lowP likelihood combination, 
we have performed an analysis including high-f polarization. In 
particular, we have considered the case of Planck TT, TE, EE+ 
lowP. The result is: B ] \ 1pc < 0.86 nG at the 95 % CL. Due to 
the nature of the effect of PMFs on the thermal history of the 
Universe and its impact on the CMB angular power spectra, the 
polarization data on small angular scales tighten the constraints 
of this analysis. 

3. Magnetically-induced non-Gaussianities 

The CMB anisotropies induced by PMFs are non-Gaussian. This 
is because magnetic forcing (as described by the magnetic en¬ 
ergy momentum tensor) is quadratic in the magnetic fields and 
therefore the resulting fluctuations are non-Gaussian even for 
Gaussian fields 10 * (Brown & Crittenden, 2005). There are al¬ 
ready published theoretical studies of the passive-mode bispec¬ 
tra (Trivedi et al., 2010; Shiraishi et al., 2011, 2012; Shiraishi, 
2013), as well as studies of the compensated-mode bispectra 

9 As explained by Chluba et al. (2015), the effect of PMFs on the ion¬ 
ization history was previously overestimated using an approach similar 
to the one of RECFAST. The reason is that in RECFAST the photoion¬ 
ization rates are evaluated assuming temperature T = T e , although from 
a physical point of view the radiation temperature should be used. This 
reduces the effect on the ionization history by up to one order of mag¬ 
nitude and a consistent treatment is implemented both in RECFAST++ 
and CosmoRec (Chluba & Thomas, 2011). 

10 This peculiarity is in common with topological defects (see 

Figueroa et al., 2010) 


(Seshadri & Subramanian, 2009; Caprini et al., 2009; Cai et al., 
2010; Shiraishi et al., 2010; Kahniashvili & Lavrelashvili, 2010) 
and of trispectra (Trivedi et al., 2012, 2014). This illustrates that 
it is possible to use CMB non-Gaussianities to constrain the 
PMF amplitude for different generation mechanisms. Several 
non-Gaussianity constraints have previously been used for this 
purpose (Caprini et al., 2009; Seshadri & Subramanian, 2009; 
Trivedi et al., 2010; Shiraishi et al., 2012; Trivedi et al., 2012). 
The non-Gaussianity constraints on PMFs are complementary 
to those derived from the angular power spectra. In this sec¬ 
tion we present three different methods for constraining PMFs 
using non-Gaussianity measurements, all involving the first of 
the higher-order stastical moments, the bispectrum. The methods 
can be applied to either the passive or the compensated modes. 

3.1. Magnetically-induced passive-tensor bispectrum 

The goal of this subsection is to derive an observational limit 
on the PMF strength from the passive bispectrum. The dominant 
contribution to the passive bispectrum is the large-scale tensor 
mode, while the scalar mode contributes subdominantly to the 
small scales. According to Shiraishi et al. (2012) and Shiraishi 
(2013), the signal-to-noise ratio (integrated over i) is expected 
to be almost saturated beyond C = 500 in estimates based on 
the temperature bispectrum. Including higher multipoles would 
therefore not bring significant improvements. Thus, we take into 
account the tensor-mode contribution for t < 500 in the fol¬ 
lowing. Here we concentrate on the almost scale-invariant case, 
n B = -2.9. In this case, the passive-tensor bispectrum is ampli¬ 
fied in the squeezed-limit configuration with l\ <s G ~ Gj, as 
a consequence of the local-type structure of the non-Gaussian 
gravitational waves induced by the PMF, given by 

h u (k) « -1.8 ln ^ TvlTB) M ijkl (k) f -^B k (p)B,(k-p). (20) 

47rp y , 0 J (2 ny 

Here r v and t b are defined in the same way as in the previous 
section, while p y o is the present photon energy density. The pro¬ 
jection tensor Mjjkitk ) is given by the products of the spin-+ 
transverse-traceless tensors as Mjjki(k) = Zj=± e^.\k) e^*(k), 
normalized as Ad,/,-/ = 4. This projection induces a tangled 
angular dependence on k , in the primordial gravitational wave 
bispectrum and the resultant CMB bispectrum is given by a 
non-factorizable combination of Gmodes (Shiraishi et al., 2011, 
2012; Shiraishi, 2012). The resultant CMB temperature and E- 
mode bispectra are almost uncorrelated with the usual scalar¬ 
mode bispectra because of their different CMB transfer func¬ 
tions. To derive constraints, we introduce an amplitude param¬ 
eter proportional to the amplitude of the magnetically-induced 
bispectrum. 


/ B 1 Mpc \ 6 

In (tJtb) 

\ 3nG / 

ln(10 17 ) 


where B\ V i pc and t b are treated as free parameters. The normal¬ 
ization factors in the last equation are chosen to be comparable 
to current upper bounds on B\ Mpc for a PMF created at the GUT 
epoch, i.e., t v /t b = 10 17 . It can be seen that the magnetically- 
induced bispectrum, which is proportional to (B i Mpc) 6 - has a log¬ 
arithmic dependence on t b . An analysis of the constraints from 
WMAP data is presented by Shiraishi & Sekiguchi (2014), yield¬ 
ing A°f a = -1.5 ± 1.4 (68 % CL). 

In order to constrain the non-factorizable magnetically- 
induced bispectrum, we use an optimal estimator derived within 
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the so-called separable modal methodology (see Fergusson et al. 
2010, Fergusson et al. 2012, Shiraishi et al. 2014, and Shiraishi 
et al. 2015 for auto-bispectra and Fergusson 2014, and Liguori 
et al. 2015 for cross-bispectra), where the theoretical bispec¬ 
trum templates are decomposed in finite subsets of the separable 
eigenbasis. Thus, the bispectrum estimator remains factorizable 
like in the usual KSW approach (Komatsu et al., 2005). Our ten¬ 
sor bispectrum template can be reconstructed well in this modal 
decomposition with about 400 eigenvectors composed of poly¬ 
nomials and a few special functions modelling the CMB temper¬ 
ature transfer functions. 


Table 4. Planck constraints on the amplitude of the non-helical 
magnetic field component, fiiMpc [nG], from the SMICA, NILC, 
SEVEM, and Commander foreground-cleaned maps at 95% 

CL. 



SMICA 

NILC 

SEVEM 

Commander 

^lMpc/flG . . 

<4.5 

< 4.9 

< 5.0 

<5.0 


pansion coefficient cl (Shiraishi et al., 2013), 


From the foreground-cleaned SMICA temperature map, we 
obtain observational constraints on the amplitude of the passive- 
tensor bispectrum. The observational data and the (Gaussian) 
simulation maps used in the computation of the linear term and 
the error bars are inpainted in the same manner as for the Planck 
tensor non-Gaussianity analysis (Planck Collaboration XXIV, 
2014; Planck Collaboration XVII, 2015), after including exper¬ 
imental aspects (beam, mask, and anisotropic noise). Our final 
result is A° ata = -1.6 ± 1.3 (T only) at 68% CL, giving no 
evidence for a signal at the 2cr level. This Planck temperature 
constraint is in good agreement with the WMAP one (Shiraishi 
& Sekiguchi, 2014). Analogous results have been derived for the 
combination of T- and A-modes and the A-only case, but since 
these results are still preliminary, we use the '/’-only mode for 
the present analysis. 

Assuming that the bispectrum is generated by PMFs, its am¬ 
plitude is given by Eq. (21). The amplitude of the bispectrum 
depends on the amplitude of the fields to the sixth power and 
on the logarithm of the ratio t v /t b , which is greater than unity. 
Therefore we have A ^ AG > 0. The result obtained in this anal¬ 
ysis therefore leads to an upper bound on the strength of GUT 
generated PMFs (ln(r v /r B ) = 10 17 ) of B\ M pc < 2.8 nG (95% 
CL). 

3.2. Magnetically-induced anisotropic passive scalar 
bispectrum 


«X*iK(* 2 K(* 3 )> = (2tt) 3 5 (3) (* 1 + k 2 + k 3 ) 

Y c L (P L (k\ ■ k 2 ) P((ki) Pc(k 2 ) + 2 perm.) , 


where c o is related to the local-form /nl as Co = 6/5 and Pi 
is the Lth order Legendre polynomial. 

If the magnetic field is generated at the GUT scale with a 
nearly scale-invariant spectrum, the Legendre coefficients are re¬ 
lated to the field amplitude via 


co 


-2 x 10“ 


^1 Mpc 

nG 


ci * -0.9 


li 1 Mpc { £>\ Mpc ^ 


nG ) \ nG 


c 2 ~ -2.8 x 10" 


B 


1 Mpc 

nG 


(23) 

(24) 

(25) 


where B\ m pc and Si Mpc are the amplitudes of the non-helical 
and helical magnetic field components (smoothed on a scale of 
1 Mpc), respectively. Estimating c 2 allows us to constrain B\ yi pc , 
but estimating ci does not lead to a useful constraint due to its 
dependence on the helical component of the PMF, which is not 
considered in this analysis (cf. Eq. 24). By the central limit the¬ 
orem, the estimated value of c 2 follows a Gaussian distribution. 
Therefore, the log-likelihood is given by 


Anisotropic stress from magnetic fields leads to curvature per- x _2 

turbations on super-horizon scales according to (Shaw & Lewis, In P{'(: 2 (d)\B\ m pc ) « — -—- - lncr, 

2010 ) 2cr 

where 


(26) 


& * 0.9 In 


1 -) 


1 


\ t b ) 4-npyfl 


§^B,(k')B j( k-k'). 


Shiraishi (2012) showed that the three-point correlation of the 
curvature perturbation £* sourced by magnetic fields is 

{((k\ ) ({k 2 ) ((kj)) oc 

PbOtP)P B ( ki)P B (k 2 )^pl 2 + /4s + /4i - \ -/h2/k!3P3ij 
~Pn(k») Pnik]) Pniki) |//23/(3i - 

+ cyclic permutations, (22) 


where p a b = k a ■ k B and k* denotes the pivot wavenumber. We can 
investigate this primordial non-Gaussianity by estimating the ex- 


x = 


c 2 (d) + 2.8 x 10 


Mpc \ 6 

”^G“/, 


with c 2 (d) being the estimated value from the data d, and <x 2 cor¬ 
responding to the variance of its estimation, which includes cos¬ 
mic variance and noise variance. In Eq. (26), we have dropped 
an irrelevant constant term. We estimate c 2 from SMICA, NILC, 
SEVEM, and Commander foreground-cleaned maps. The vari¬ 
ance is estimated from the realistic Planck simulations for each 
foreground-cleaning method (Planck Collaboration IX, 2015; 
Planck Collaboration X, 2015). To determine B\ m pc and its con¬ 
fidence region, we use the CosmoMC package (Lewis & Bridle, 
2002) as a generic sampler and obtain the posterior probability 
of B\ Mpc, given the likelihood. This analysis yields upper bounds 
on the amplitude of the non-helical magnetic field component, 
B i m P c, which are presented in Table 4. Using the Planck data 
constraint on the local-form /nl (Planck Collaboration XVII, 
2015) and co = 6/5 /^j aal , we impose an additional constraint 
of B i Mpc < 5.5 nG at 95 % CL, which is weaker than the C 2 con¬ 
straints. 
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The constraint from the passive-scalar bispectrum with 
SMICA maps is .SiMpc < 4.5 nG (95% CL). Thus, the ad¬ 
dition of the polarized bispectrum leads to an improved con¬ 
straint with respect to the previous limit, B\ yi pc < 5.2 nG, 
from the Planck 2013 analysis (Shiraishi et al., 2013; Planck 
Collaboration XXIV, 2014). 


3.3. Magnetically-induced compensated-scalar 
bispectrum 

Now we derive the magnetically-induced scalar bispectrum on 
large and intermediate angular scales using a semi-analytical 
method. We compute an effective / NL based on the comparison 
between the bispectrum and the power spectrum and derive the 
constraints on the amplitude of the PMF using Planck measure¬ 
ments. 

We derive the magnetically-induced scalar bispectrum on 
large angular scales for compensated initial conditions, basing 
our analysis on the treatment presented by Caprini et al. (2009). 
For simplicity we redefine the parameter describing the ampli¬ 
tude of the PMF in this section. Instead of using the smoothed 
amplitude we directly use the root mean square value of the field. 
This quantity is finite thanks to the sharp cut-off inserted in the 
PMF power spectrum to model the small-scale suppression of 
the field. The mean square of the field is then defined as 


(B\k)) = 


A b k’^ 

2n 2 n B + 3 


(27) 


The magnetically-induced bispectrum on large angular scales 
depends on the temperature anisotropy on large angular scales 
and therefore on the Sachs Wolfe signal induced by the PMFs 
(Caprini et al., 2009; Bonvin & Caprini, 2010). We use the 
expression derived by Paoletti et al. (2009) and Caprini et al. 
(2009), 

©f (r 0 , k) 

e 2f+l * orO B (k) j f (k(To - r dec )), (28) 


The characteristic feature of the magnetically-induced bis¬ 
pectrum, generated by compensated modes, is that, contrary to 
what often happens for inflationary non-Gaussianities, it is not 
possible to identify an a priori dominant geometric configura¬ 
tion. It is therefore necessary to analyse the bispectrum inde¬ 
pendently of the geometric configuration. Caprini et al. (2009) 
derive an approximate expression for the three-point correlation 
function of magnetic energy density, which is independent of 
the geometric configuration and is tested against the analytic re¬ 
sults for the flattened case. Using this expression, we derive the 
magnetically-induced bispectrum and specify a geometric con¬ 
figuration only after integrating the magnetic energy density bis¬ 
pectrum in k-space. We show the results for the case correspond¬ 
ing to a local /nl- The magnetically-induced bispectrum for the 
nearly scale-invariant case, hb = -2.9, is given by 
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(31) 


For the specific case n B — -2 it is 


j'eff 

/nl 


57T u a 3 (B 2 y It, 


= 5.19 


P rel 

(B 2 ) 


(10 9 G) 2 


log(lp/4iax) 

log(4iax/4iin) 


(32) 


where 0^ o> is the temperature anisotropy, il B (k) = (B 2 (k)) /p Te \, 
and To and Td ec are the conformal time at present and at decou¬ 
pling, respectively. For simplicity we have used an approximated 
expression for the initial conditions instead of the exact one, 
which involves also the Lorentz force (Caprini et al., 2016). We 
therefore introduce a correction factor, a = 0.5, which numeri¬ 
cally includes all the contributions that we do not consider in the 
expression. The above relation holds for the compensated mode 
initial conditions. The CMB bispectrum on large scales can then 
be written as 

(47t) 3 (-0 fl+4+f3 

{a (imi a tim2 a {3m ,) - (2 ^ + i)( 2 £, + l)(2/ 3 + 1) X 
r d\ d 3 q d 3 p 

J (27r) 9 

KmSbrUiqiYUiP) 

x<0j(To, k) ®f(ro, q) (to, p)) . (29) 

Substituting Eq. (28), we note that the bispectrum depends on 
the 3-point correlation function of magnetic energy density, 

, n\ r \ r w 1 f d 3 kd 3 qd 3 p 

(Bi(k) B/(k - k) Bj(q)Bj(q - q) B,(p) B,(p - p)). 


where we assume Jt = 18.98 x 10 9 as the amplitude of the pri¬ 
mordial gravitational potential power spectrum ( t 2 Cc = Jt/n). 
In all numerical estimates we have taken Hu = k n r (l = 3000, 
/ max = 750, and ( nlm = 10. The Planck limit on local / NL 
is given in Planck Collaboration XVII (2015). We use the re¬ 
sult of the SMICA KSW T+E ISW-lensing-subtracted analy¬ 
sis, namely / NL < 5.8 at 68 % CL. This limit on / NL translates 
into constraints on the PMF amplitude of \/(B 2 ') < 1.8 nG for 

n B = -2.9, VCS 3 ) < L0 nG for n B = -2, and ^|{B 2 ) < 1.7 nG 
for n B = 2. The constraints on the smoothed amplitude of the 
field are fiiMpc < 3nG for n B = -2.9, Bim P c < 0.07 nG for 
n B = -2, and B\ vi pc < 0.04 nG for n B = 2. These results show 
how the constraints are competitive with, in addition to being 
complementary to, the ones given by the CMB angular power 
spectrum. 


4. Faraday rotation 

4. 1. Constraints on PMFs from the Faraday rotation 
power spectrum 

The presence of a PMF at the last scattering surface induces a ro¬ 
tation of the polarization plane of the CMB photons (Kosowsky 
& Loeb, 1996). This effect is known as Faraday rotation (here¬ 
after FR). The Faraday depth <1> is proportional to the integral 
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along the line of sight of the magnetic field component along 
this direction, By, and the thermal electron density, « e , i.e., 


ch = K 


f 


n e (x, h) ft) dx. 


(33) 


The constant is K = 0.81 radm^pc -1 cm 3 /iG _1 = 2.6 x 
10~ 26 radnG 1 . The unit of the Faraday depth is rad m 2 . In 
this equation, h is the unit vector in the line-of-sight direction. 

Here we analyse the impact of FR on the polarized CMB 
power spectra. We assume that the magnetic field is generated 
at some pre-decoupling epoch. We do not consider the genera¬ 
tion mechanism itself, but only investigate the observable effects 
caused at recombination. Several works have derived the modi¬ 
fication of the Boltzmann equation for the Stokes parameters in 
the presence of a homogeneous PMF (see, e.g., Scoccola et ah, 

2004) and for a stochastic distribution (see, e.g., Kosowsky et ah, 

2005) . Here we explore the case of a stochastic distribution. 

As discussed above, a PMF may induce scalar, vector, and 
tensor perturbations. At recombination, FR mixes the signatures 
of different perturbations. Some previous attempts to constrain 
FR signatures in the WMAP power spectrum are presented by 
Kahniashvili et ah (2009) and Pogosian et ah (2011). Both as¬ 
sume that /-’-modes are converted into /(-modes via FR. They 
obtain magnetic field strength limits of B\ m pc $ lOOnG and 
suggest an almost scale-invariant spectmm, i.e., «b ~ -2.9 for a 
power-law distribution. 

As in previous works, we assume a PMF distribution de¬ 
scribed by a power law as in Eq. (2). Note that any helical part of 
the field does not contribute to the Faraday rotation (Campanelli 
et ah, 2004). The generation of magnetically-induced /(-modes 
through FR of E-modes is described by (Kosowsky et ah, 2005) 


s-i/BB 




tih 


(2Ei + \)(2h + 1) 

4;r(2 £ + 1) 


Nlm,c u { 2 fcl E ci (c f %) 2 


(34) 


and the rotation of primordial /(-modes into magnetically- 
induced E-modes by 


sirEE 


2 v (2R + l)(2l 2 + 1) 


x 


NlWM,h?CfCl (C f %) 2 


(35) 


where N { = (2(£ - 2)!/(E + 2)!) 1/2 is a normalization factor, 
K{£,{\,{ 2 ) = -1/2 (E 2 + E 2 + E;-2EiE 2 -2E,E + 2E[ -2E 2 -2E) 
with E = /’(/’+ 1), Ei = C\{£\ + 1), L 2 = C 2 (f 2 + 1), and C® 0fi0 is a 
Clebsch-Gordan coefficient. The power spectrum of the rotation 
angle is related to the one of the Faraday depth through 

C a e = Vo 4 C®, (36) 


where vq is the observed frequency, and 



9({(+\) ^1 Mpc 

(47r) 3 e 2 r(« B + 3/2) 



dxx J,B j 2 f (x). 


(37) 


Here, xd = r k n , where r is the conformal time and k\> is given 
by Eq. (3). 

In Eqs. (34) and (35), C EE and C BB are the primordial power 
spectra, whereas C' ( EE and C' BB are the ones including the effect 
of Faraday rotation, i.e., the observed ones. We use the observed 
E-mode spectrum C' EE at 70 GHz as a proxy for the primordial 


o 
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Fig. 12. Probability contours of PMF strength vs. spectral index 
of the PMF power spectrum as constrained by the 70 GHz obser¬ 
vations. 


one, C EE , in Eq. (34) to calculate predicted /(-mode spectra to 
be compared with the observed one. 

The 70 GHz observations give /([ Mpc < (1040, 1380) nG 
(68%, 95% CL). The reduced^ 2 is 1.35. The magnetic field 
spectral index for the stochastic distribution remains uncon¬ 
strained. In Fig. 12, we show the probability contours derived 
from the 70 GHz data for the magnetic field strength /(1 Mpc and 
the spectral index of the PMF power spectrum, The upper 
bounds obtained from the FR analysis are very high compared to 
the other methods. They are also slightly weaker than those pre¬ 
viously obtained from FR analyses with WMAP-5 (Kahniashvili 
et ah, 2009), WMAP-7 (Pogosian et ah, 2011), and WMAP-9 
(Ruiz-Granados et ah, 2015) data. 

4.2. Robustness of the results in the presence of 
foregrounds 

The magnetic fields of the Milky Way contribute to the net 
Faraday rotation. Although the precise geometry of these mag¬ 
netic fields remains uncertain (see, e.g., Ruiz-Granados et ah 
2010 and lansson & Farrar 2012), the Galactic Faraday depth 
could be a foreground for the primordial Faraday depth, at least 
on large scales. 

To quantify the impact of the Galactic FR on the detection 
of primordial magnetic fields, we use Galactic observations of 
polarized synchrotron emission at 1.4 GHz and 23 GHz and the 
synthesized all-sky Faraday rotation map derived from extra- 
galactic radio source emission provided by Oppermann et ah 
(2014). In addition, we use simulations of the Galactic Faraday 
rotation obtained by using an axisymmetric Galactic magnetic 
field model for the halo field described by Ruiz-Granados et ah 
(2010). Maps of Stokes Q and U are provided by Wolleben et ah 
(2006) at 1.4GHz and by Bennett et ah (2013) at 23GHz. Both 
frequencies are dominated by polarized synchrotron emission 
from within the Milky Way and are used to obtain the Galactic 
Faraday depth (see Ruiz-Granados & Florido, 2015, for details). 
For computing the power spectrum of the FR coming from sim¬ 
ulations and observations at 1.4 and 23 GHz, we use the polar¬ 
ization proccessing mask provided by WMAP-9 11 . 


11 http://lambda.gsfc.nasa.gov/product/map/dr5/ 
m_products.cfm 
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Fig. 13. Power spectra of the primordial Faraday depth <1* for 
magnetic field strengths of 100 (solid line), 10 (dotted line), and 
1 nG (dashed line), the Galactic Faraday depth from the all-sky 
Faraday map of Oppermann et al. (2014) (blue dashed line), 
the Galactic Faraday depth from Oppermann et al. (2014) at 
45° latitude (green dashed line), the Galactic Faraday depth de¬ 
rived from Galactic emission at 1.4 and 23 GHz (red circles), and 
the Galactic Faraday depth for a Galactic magnetic field model 
(black crosses). 


In Fig. 13, we show the power spectra of Galactic <I> de¬ 
rived from polarized measurements at 1.4 GHz and 23 GHz, 
simulations, and for the all-sky Faraday rotation map provided 
by Oppermann et al. (2014) 12 . The fluctuations in the Galactic 
Faraday sky are not isotropic. Therefore, their statistics are not 
completely described by a power spectrum. Oppermann et al. 
(2014) model the Galactic Faraday depth as the product of an 
isotropic Gaussian random field and a latitude-dependent func¬ 
tion. In Fig. 13, we show two angular power spectra derived from 
the results of Oppermann et al. (2014). For the first one, we 
generate Gaussian realizations from their angular power spec¬ 
trum, multiply them with their latitude profile, and pass them 
through the anafast routine of HEALPix. Averaging the re¬ 
sult over 1000 realizations gives the blue dashed line in Fig. 13. 
For comparison, we also show the angular power spectrum of 
Oppermann et al. (2014) multiplied with the square of their pro¬ 
file function at a latitude of \b\ = 45°, which gives the strength 
of the foreground Faraday rotation at a typical latitude used for 
CMB analysis. We plot also the power spectrum of the primor¬ 
dial Faraday rotation for PMFs of strength 1, 10, and 100 nG, 
respectively, and a spectral index of n B = -2.9. 

Our main conclusion, as indicated by Fig. 13, is that obser¬ 
vations of the Galactic Faraday depth show that primordial mag¬ 
netic field strengths lower than lOnG would require a detailed 
knowledge of the Galactic FR, at least at multipoles lower than 
C ~ 50, where the foreground rotation would be dominant. The 
exact detection limit for PMFs achievable through a study of 
Faraday rotation depends on Galactic latitude and on the extent 
to which the correction for foreground rotation is possible. 

The upper limit found from the FR analysis is B\ m P c < 
1380nG and no restriction can be obtained for the spectral in¬ 
dex n B . In any case, our FR constraint on PMFs is well above 
the expected contamination level by Galactic FR, and therefore 

12 http://www.mpa-garching.mpg.de/ift/faraday 


this contamination is currently not an issue. This is consistent 
with the prediction by De et al. (2013) that the Galactic Faraday 
depth would not be measurable with Planck. 


5. Constraints on PMFs from Alfven waves 


Here we investigate the signature statistical anisotropy induced 
by Alfven waves, which delivers yet another constraint on PMFs. 
In Planck Collaboration XXIV (2014) we constrained the Alfven 
waves in the early Universe, where some arbitrary origin (includ¬ 
ing stochastic PMFs) was assumed for primordial vector pertur¬ 
bations. Given no evidence for Alfven waves from that analysis, 
we now consider stochastic PMFs as the source of primordial 
vector perturbations and constrain an average background mag¬ 
netic field and the energy density of stochastic PMFs 13 . PMFs 
may produce Alfven waves in the early Universe, which leave 
observable imprints on the CMB via the Doppler and integrated 
Sachs-Wolfe effects. Durrer et al. (1998) show that Alfven waves 
in the early Universe generate a fractional CMB anisotropy 

AT 

—(«, *) ~ n • kl(k, r last ) = n ■ Sl 0 (k) v A k Ti ast h 0 ■ k , (38) 

i o 

where k denotes a Fourier mode vector, h a sky direction, ho 
the unit vector in the direction of the homogeneous background 
magnetic field B, and Tq is taken again as 2.7255 K (Fixsen, 
2009). Here £l(k,r i ast ) and Hoik) denote the gauge invariant lin¬ 
ear combination of vector perturbations at last scattering and at 
an initial time, respectively. In this analysis, we assume a non¬ 
helical stochastic PMF, B, to be the sole source of initial vec¬ 
tor fluctuations, £io = v A /B B (Durrer et al., 1998; Kahniashvili 
et al., 2008). The Alfven wave velocity, v A , is given by (Durrer 
et al., 1998) 


v A 


B 

2 y/n(p r + p r ) 


» 2.2 x 10 5 m s _1 


B 

TnG ’ 


(39) 


where p r and p r are the co-moving density and pressure of the 
photons. 

Kahniashvili et al. (2008) show that Alfven waves in the 
early Universe produce correlations between harmonic modes 
separated by AC = 0,±2, and Am — 0, ±1,±2. We give the ex¬ 
plicit form of the correlations in Appendix B. Investigating these 
imprints, we impose a constraint on the Alfven waves in the 
early Universe. In the weak Alfven wave limit, the CMB data 
log-likelihood £. can be expanded as 


r n d£ - 
-L ~ -L\h=o + o. 


dli 


h=0 


, \d 2 L 
1 + 2 dh 2 


-0(h 3 ), 


h=0 


where h — B 2 Mpc u 3 / B 2 . The first term on the right hand side 
is simply equal to the likelihood of the standard cosmological 
model and the first and second derivatives of the likelihood are 
obtained by 

df d 2 f 

-=■ = Pi - (<H), — =■ = -CW 2 > + , (40) 


where (...) denotes the ensemble average of signal and noise. 
The Hessian TT is given by 




(41) 


13 For this analysis we follow the approximation of considering the 
stochastic background of PMFs as split into an average background 
field, which emulates the average effect of the stochastic fields, and the 
stochastic fields at the perturbative level. 
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Table 5. Planck constraints on the Alfven wave amplitude 
BjMpc/5 2 v a- The Alfven wave velocity ua is normalized to the 
speed of light. 


Confidence Level 

68 % 

95% 

fimpc IB 2 < . 

< 3.4 x 10- 7 

< 1.7 x 10~ 5 


where a is the vector consisting of the spherical harmonic coef¬ 
ficients, a ( m , of the CMB anisotropy data, and C is their covari¬ 
ance matrix. 

In our analysis, we consider the foreground-cleaned SMICA 
map, where we apply the common mask (Planck Collaboration 
X, 2015). We assume the fiducial Planck cosmological model 
and use realistic Planck simulations to estimate the ensemble 
average values for signal and noise, as required in Eq. (40). The 
quantity C [ a, required in Eq. (41), is determined via the mes¬ 
senger field method (Eisner & Wandelt, 2013). Some of the pa¬ 
rameters (hb,6b,<Pb) influence the signature correlation nonlin- 
early. Due to these nonlinear parameters, we use the CosmoMC 
package (Lewis & Bridle, 2002) as a generic sampler for the log- 
likelihood in Eq. (40) and obtain the posterior probability for the 
Alfven wave parameters \B 2 lMpc /B 2 v A A , n B , 0 B , <Pb)- As discussed 
previously, we assume the initial vector fluctuations to be en¬ 
tirely sourced by a non-helical stochastic PMF, B. In Table 5, we 
show upper bounds on this combination of parameters at 68 % 
and 95 % CL, after marginalizing over the spectral index n b and 
the direction 0 B , <Pb- 

Other theoretical models with correlations across multipoles 
with AC = ±1, ±2 are investigated in Planck Collaboration XVI 
(2015) and Planck Collaboration XX (2015). The Planck data 
show no evidence in favour of these models. 


6. Conclusions 

6.1. Methodology 

In this paper, we have presented constraints on a stochastic back¬ 
ground of primordial magnetic fields using Planck data. PMFs 
may have left different types of imprints on the CMB. This is 
why the CMB can be regarded as one of the best laboratories 
for investigating and constraining PMFs. The richness of Planck 
data, which provide several different probes based on the statis¬ 
tics of CMB anisotropies, allows us to constrain PMFs using 
different methods, deriving constraints that are complementary 
to each other. Aiming for a broad perspective on PMFs, we have 
taken advantage of these different possibilites offered by Planck. 
In particular, we have exploited the impact of PMFs on the CMB 
anisotropy angular power spectra in temperature and polariza¬ 
tion, both through the magnetically-induced modes and their ef¬ 
fect on the ionization history of the Universe. In addition we 
have considered the Faraday rotation induced by them on the 
CMB polarization. Beyond the two-point statistics probed by 
the angular power spectrum, we have investigated higher-order 
statistical moments of the CMB imprinted by PMFs. In particu¬ 
lar, we have analysed the CMB bispectrum for the presence of 
magnetically-induced passive and compensated modes. Finally, 
we have considered the correlations between harmonic modes 
induced by PMFs. 


6.2. Constraints for non-helical fields from the angular 
power spectra 

The constraints based on the CMB angular power spectra have 
been derived using the Planck likelihood. The general analy¬ 
sis, which considers only the contributions from compensated 
modes, provides the constraint B\ m pc < 4.4 nG at 95 % CL, with 
positive spectral indices constrained to lower amplitudes of the 
fields, which we have shown to be robust under the inclusion 
of high-/: polarization data. In fact, the impact of PMFs on the 
EE and TE polarization spectra is negligible compared to the 
dominant contribution given by the magnetically-induced vec¬ 
tor modes on the TT spectra on small angular scales. The in¬ 
clusion of the passive tensor contribution does not improve the 
constraints on the amplitude but affects the posterior distribution 
for the PMF spectral index. The results confirm that CMB data 
constrain positive spectral indices to smaller amplitudes. The in¬ 
clusion of the contribution of passive modes, thanks to their sen¬ 
sitivity to the slope of the PMF spectrum, strongly disfavours 
nearly scale-invariant magnetic spectra for amplitudes greater 
than very few nanogauss (see Fig. 8). This is due to the large con¬ 
tribution that magnetically-induced passive tensor modes with 
nearly scale-invariant spectra give to the CMB anisotropies on 
large angular scales. 

We have also performed an analysis using both the Planck 
2015 data and the BICEP2/ Keck-Planck cross-correlation. The 
results are fully compatible with the analysis based on Planck 
data alone, with only a slightly higher upper limit on the PMF 
amplitude. 

Our likelihood analysis for PMF properties is sensitive to 
CMB foreground residuals in the data, since these contribute to 
small angular scales in the CMB spectra as well. In our analy¬ 
sis we have used the foreground residual treatment provided by 
the Planck likelihood. Since some of the foreground models as¬ 
sume an angular power spectrum with a shape similar to the one 
given by magnetically-induced perturbations, we have investi¬ 
gated the issue of possible degeneracies between PMFs and fore¬ 
ground residuals. In particular, we have noticed that there is a de¬ 
generacy with the Poissonian terms modelling unresolved point 
sources at 100, 143, 143x217 and 217 GHz. Fortunately, other 
foreground contributions do not show any degeneracy thanks 
to their different spectral shapes. We have further tested how 
severely foreground residuals might affect the obtained con¬ 
straints on PMF amplitudes. Since PMFs that are consistent with 
the data do not significantly affect other cosmological parame¬ 
ters, we have fixed the amplitudes of foreground residuals to the 
values obtained under the Planck ACDM cosmology. The limit 
in this case, B\ m rc < 3.0 nG at 95 % CL, is slightly tighter than 
the limit obtained without fixing the foreground residual ampli¬ 
tudes, B i m P c < 4.4 nG at 95% CL. This test has no statistical 
significance in the constraints on PMFs but illustrates the im¬ 
pact foreground residuals can have on constraints on the PMF 
amplitude. 

Our new constraints are compatible with previous constraints 
from other experiments and from the previous Planck release 
(Planck Collaboration XVI, 2014). The slightly higher upper 
limits with respect to the 2013 Planck release are due to changes 
in the 2015 data. In particular, the changed calibration and 
slightly different slope of the power spectrum of cosmological 
perturbations allow for stronger PMFs, with possible contribu¬ 
tions from the different foreground residual treatment. 
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6.3. Constraints on maximally helical PMFs 

We also constrain maximally helical PMFs. We restrict our anal¬ 
ysis to the maximally helical case because of the absence of 
TB and EB information in the Planck 2015 high-/ - likelihood. 
Maximal helicity decreases the amplitude of the magnetically 
generated CMB fluctuations and, as a consequence, we obtain 
B \Mpc <5.6 nG at 95 % CL in this case. 

6.4. Selected scenarios 

We have further investigated two specific PMF models of inter¬ 
est: causally generated fields with a spectral index of hb = 2 
and fields with an almost scale-invariant power spectrum with 
n B — -2.9. The constraints for these extreme cases are B"^ 2 c < 
0.011 nG and /('"Q; 2<J < 2.1 nG at 95 % CL, respectively. 

The impact of PMFs on the ionization history of the Universe 
directly affects the CMB temperature and polarization power 
spectra. In particular, we have also considered the main dissi¬ 
pative effects operating during and after recombination, namely 
ambipolar diffusion and energy cascading in MHD turbulence 
in the prediction for the CMB spectra in temperature and polar¬ 
ization. These modify the primary CMB power spectra in ad¬ 
dition to the gravitational contributions of the magnetic modes. 
For the nearly scale-invariant case we have obtained the con¬ 
straint B"m pc 2 9 < 0-9 nG at 95 % CL. This limit is tighter than 
when neglecting the effect on the ionization history. However, 
uncertainties related to the modelling of the heating mechanism 
(see discussion by Chluba et al., 2015) suggests that further in¬ 
vestigation of this promising avenue is needed. 

6.5. Non-Gaussianity-based constraints 

For the non-Gaussianity analyses we have focused on the passive 
modes with a nearly scale-invariant power spectrum, n B — -2.9, 
and the compensated scalar modes. These are the dominant con¬ 
tributions on large angular scales, where the non-Gaussianity 
analyses are performed. 

In our first CMB non-Gaussianity analysis, we have consid¬ 
ered passive tensor modes for PMFs with nearly scale-invariant 
spectra. These contribute predominantly to the CMB fluctua¬ 
tions on large angular scales. For this case, we have calculated 
the resulting CMB bispectrum and compared it with the obser¬ 
vational limit. We have used a bimodal decomposition to esti¬ 
mate the amplitude of the Planck bispectrum in the squeezed 
configuration, in which the observational limit on the amplitude 
of the bispectrum can be translated into a constraint on the am¬ 
plitude of PMFs. Using the temperature bispectrum we obtain 
B'lVf < 2.8 nG for fields that were generated at the GUT 
phase transition. 

For our second non-Gaussianity analysis we have used a 
different approach to the magnetically-induced bispectrum. We 
have considered the passive contributions by tensor and scalar 
modes for nearly scale-invariant fields, but instead of using 
the bispectrum amplitude we have used the local type of non- 
Gaussianity induced by PMFs, considering only the non-helical 
part. The local bispectrum from multipoles with X 2 ] ( n = even 
(where l n stands for the multipoles l\,G, (3 normally used to 
express the bispectrum, see for example Eq. 29) has been decom¬ 
posed with coefficients determined by the amplitude of the PMF. 
A likelihood estimation allows us to constrain the PMF am¬ 
plitude with CMB foreground-cleaned maps. We have applied 
this method to all four component separation methods used by 


Planck. The SMICA map, which is expected to contain the least 
foreground residuals (Planck Collaboration IX, 2015), gives the 
constraint B"^ 2 ' 9 < 4.5 nG, improving previous constraints de¬ 
rived from the scalar bispectrum for WMAP data. 

Our third non-Gaussianity analysis focuses on the compen¬ 
sated scalar modes. In this case, we have used an analytic esti¬ 
mate of the bispectrum on large angular scales. We have used 
an improved estimate of the source term (Caprini et al., 2016) 
with respect to previous results (Caprini et al., 2009). This ana¬ 
lytic estimate can be compared with the observed local /nl from 
Planck, giving B" B ^ 2 ' 9 < 3.0 nG. 

The results from the different non-Gaussianity analyses (al¬ 
though coming from different methods) are all consistent and at 
the level of those derived with the likelihood analysis using only 
the CMB angular power spectra. 

6.6. Constraints from Faraday rotation 

We have further considered the effects of Faraday rotation on 
the primary CMB polarization anisotropies. In this context, we 
have used the EE- and /(/(-polarization power spectra. We have 
derived the constraints on the PMF amplitude using a y 2 anal¬ 
ysis based on the LFI 70 GHz low -/ 1 (£ < 30) polarization 
power spectra. The resulting constraint is B\ M p C < 1380 nG. 
The upper limits from Faraday rotation are larger than those de¬ 
rived from magnetically-induced perturbations, thermal effects, 
and non-Gaussianity. On one hand, the Faraday rotation signal 
rapidly vanishes with increasing frequency (see Eq. 36) and thus 
is strong only for lower frequencies. On the other hand, since the 
BB spectrum is the results of the rotation of the EE spectrum, it 
has a stronger contribution on smaller angular scales. Our anal¬ 
ysis includes only the low multipoles of the 70 GHz data, where 
the signal is lower due to its spectral shape. The combination of 
data sets and type of signal results in less stringent constraints. 
However, even with the restricted subset of Planck data avail¬ 
able, the constraints are only slightly weaker than derived in 
previous analyses performed with WMAP (Kahniashvili et al., 
2009; Pogosian et al., 2011; Ruiz-Granados et al., 2015). 

To estimate the impact of Galactic Faraday rotation on the 
results, we have analysed synthetic Galaxtic Faraday maps as 
well as radio synchrotron data at 1.4 GHz and 23 GHz. We have 
accounted for the fact that the signal is not isotropic on the sky 
but depends on the latitude of the observations. We have derived 
an estimate of the power spectrum for the Faraday depth, shown 
in Fig. 13, and compared it with the predictions for different val¬ 
ues of the PMF amplitude. Our results show that the threshold 
for which the Galactic contamination may become relevant is 
around 10 nG. This amplitude is much below our current con¬ 
straints, which can therefore be considered clean from Galactic 
contamination. 

6.7. Constraints on Alfven waves 

To complete the round of different types of analyses involving 
different probes, we have investigated the correlation induced 
between different modes in harmonic space by Alfven waves 
produced by the presence of PMFs (Kahniashvili et al., 2008). 
This correlation has not been used to constrain the PMF ampli¬ 
tude directly, but the Alfven wave parameter, which is a com¬ 
bination of the stochastic background amplitude, the Alfven ve¬ 
locity, and the assumed mean background field. We have again 
used the SMICA foreground-cleaned map to derive the upper 
limit B] Mp c v 4 a /B 2 < 1.7 x 10 5 . From this constraint we de- 
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duce that the data do not show any evidence of Alfven waves, 
a conclusion that was also reached in previous analyses (Planck 
Collaboration XXIV, 2014) carried out with more generic as¬ 
sumptions on the origin of the Alfven waves. The absence of 
Alfven waves is also compatible with the results from other mod¬ 
els for the harmonic space correlations that are not related to the 
PMFs (Planck Collaboration XVI, 2015; Planck Collaboration 
XX, 2015). 

6.8. Concluding summary 

The results presented show that CMB anisotropies are one of 
the best probes for investigating the nature of PMFs. The Planck 
2015 data offer the possibility to use the PMF’s signatures either 
in the angular power spectra in temperature and polarization or 
in higher-order statistics, where both measurements can be tack¬ 
led by different methodologies. All the independent constraints 
we obtain are consistent with each other. 

The Planck 2015 data constrain the PMF amplitude at the 
nanogauss level. Different signatures are sensitive to different 
contributions and may be optimal for specific types of PMF. 
In particular, the analysis that uses the gravitational impact of 
PMFs on the CMB angular power spectra is dominated by the 
compensated vector contribution on small angular scales and 
therefore is able to constrain PMFs without any assumptions on 
their generation mechanism. On the contrary, two of the three 
analyses of non-Gaussianities are dominated by passive tensor 
modes, which can provide significant constraints only for nearly 
scale-invariant PMFs. 

The future of both classes of methods, the angular power 
spectra and the non-Gaussianities, is bright, but three avenues 
are particularly promising. The helicity of PMFs will be con¬ 
strained by TB and EB cross-correlations, which will be in¬ 
cluded in the next Planck release. The study of the PMF’s im¬ 
pact on the ionization history is expected to further improve with 
future Planck polarization data. Non-Gaussianities are a distinc¬ 
tive signature of PMFs and further studies may provide more 
and more refined predictions of the magnetically-induced pas¬ 
sive and compensated CMB bispectra and trispectra, which will 
improve the future Planck analyses. 
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Fig. A.l. Probability contours for the PMF amplitude and the 
foreground parameters for the Planck 2013 likelihood. 


Yamazaki, D. G., Ichiki, K., Kajino, T., & Mathews, G. J., Effects of a primordial 
magnetic field on low and high multipoles of the cosmic microwave back¬ 
ground. 2008, Phys. Rev. D, 77, 043005, arXiv:0801.2572 
Yamazaki, D. G., Ichiki, K., Kajino, T., & Mathews, G. J., New constraints on the 
primordial magnetic field. 2010, Phys. Rev. D, 81, 023008, arXiv:1001.2012 
Zizzo, A. & Burigana, C., On the effect of cyclotron emission on the spec¬ 
tral distortions of the cosmic microwave background. 2005. New A, 11, 1, 
arXiv:astro-ph/0505259 


Appendix A: Impact of foregrounds on PMF 
constraints from the angular power spectra 

In Fig. A.l we present the two-dimensional proability distri¬ 
butions of PMF amplitude and foregroun d parameters for the 
Planck2013 likelihood, which show only a mild degeneracy with 
the Poissonian amplitude for the 143 GHz. In Fig. A. 2, we plot 
the two-dimensional probability distributions of the PMF am¬ 
plitude and the foreground parameters (for their description see 
Planck Collaboration XI 2015), except for the Poissonian terms, 
which have been discussed in Sect. 2.4.4. These plots do not 
show any degeneracy. 
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Fig. A.2. Probability contours for the PMF amplitude and the foreground parameters. 
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Appendix B: Statistical anisotropy induced by 
Alfven waves 

It has been shown that the presence of Alfven waves in the early 
Universe leads to specific correlations of the CMB in harmonic 
space (Kahniashvili et al., 2008). The signature correlations in¬ 
duced by Alfven waves are as follows: 
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Here Q is the power spectrum in the absence of Alfven waves, 
6 b and (b B are the spherical angles of the direction of the back¬ 
ground magnetic field B, and I f J is given by 
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where k\ } is the co-moving wave number of the dissipation scale 
due to photon viscosity and given approximately by 10/(cri ast ) 
(Durrer et al., 1998). The dissipative damping effect becomes 
significant at multipoles I > 500 (Durrer et al., 1998). We note 
that the damping scale considered in this context is different 
from the one considered in Sect. 2. The damping effects con¬ 
sidered in each case are related to two different aspects and are 
specific to the two topics treated; in the study of the impact of 
PMFs on CMB anisotropies, the damping scale considered is 
due to the dissipation of the PMFs themselves as investigated by 
Subramanian & Barrow (1998a), in this section instead, we con¬ 
sider a damping scale derived from the damping of the vector 
perturbations generated by the PMFs, not the PMFs themselves. 
The latter damping scale is derived by Durrer et al. (1998). In 
Eq. (B.l), B 2 x denotes the total power of the non-helical PMF 
smoothed at the spatial scale of A. There are general cases for 
dipole and quadrupole coupling, where cosmological parame¬ 
ters are assumed to vary with position (Moss et al., 2011). If we 
allow a position-dependent parameter to be a vector and treat the 
PMF as a position-dependent parameter, the correlation induced 
by Alfven waves may be incorporated into the framework of this 
approach. 
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